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Abstract 

Consider four point particles with equal masses in the euclidean space, subject to the 
following symmetry constraint: at each instant they are symmetric with respect to the dihedral 
group D2, that is the group generated by two rotations of angle n around two orthogonal axes. 
Under a homogeneous potential of degree —a for < a < 2, this is a subproblem of the four- 
body problem, in which all orbits have zero angular momentum and the configuration space 
is three-dimensional. In this paper we study the flow in McGehee coordinates on the collision 
manifold, and discuss the qualitative behavior of orbits which reach or come close to a total 
collision. 

MSC Subject Class: Primary 70F10; Secondary 37C80. Keywords: Dihedral 4-body problem, 
McGehee coordinates, heteroclinics. 

1 Introduction 

The goal of this paper is to investigate the qualitative behavior of solutions of the dihedral sym- 
metric 4-body problem in space under the action of a homogeneous potential of degree —a. For the 
Newtonian potential this problem is a kind of generalization of Devaney planar isosceles three body 
problem |Dev801 IDev81j . following Moeckel's approach to the study of the three body problem 
in space |Moe81[ IMoe83j . The dihedral four-body problem is a subproblem of the full four-body 
problem which reduces to a three-dimensional configuration space. Briefly, one takes n = 4 equal 
masses whose initial position and velocity are symmetric with respect to the 2-dihedral group 
(which is isomorphic to the Klein group Z 2 x Z 2 ) D 2 C SO(3). So the masses form a (possibly 
planar, degenerate and non-regular) tetrahedron in space. Because of the symmetry of the prob- 
lem, the masses will remain in such a configuration for all time. Hence we have a system with 
only three degrees of freedom. 

We will use McGehee coordinates |McG74j in order to study the dynamics of the dihedral four 
body problem for a general homogeneous potential of degree —a and with a slight change: we 
consider McGehee coordinates not only for studying the behavior of solutions passing close to a 
total collision, but also for parabolic orbits connecting central configurations, projecting the full 
phase space to a codimension 1 subspace. We replace the singularity due to total collapse with an 
invariant immersed manifold in the full phase space usually called total collision manifold which 
is the immersion of the parabolic manifold of the projected phase space. By studying this flow we 
are able to establish some global results on the behavior of solutions. We discuss the qualitative 
behavior of orbits which reach or come close to the total collision and the behavior of orbits which 
start from total collapse, which implies chaotic behavior. Such of problem includes a some of other 
subsystems with one or two degrees of freedom studied in the past years. We also observe that 
the behavior of the bodies is the same for all values of the parameters a £ (0, 2) and not only for 
the Newtonian case a = 1. 

The literature on this problem is quite broad and for this we quote just few papers in which 
similar studies were carried over. In particular some properties of the behavior of the flow on the 
total collision manifold has been established for the planar case in |SL82j . Moreover, in the spatial 
case in a series of papers Delgado and Vidal in |DV99|. IVid99j studied the so-called tetrahedral 



four body problem without and with rotation. We observe that the potential in this case (up to 
a multiplicative constant depending on the normalization of the masses) coincides with a subcase 
of the dihedral four body problem. 

The main results are presented in the first part of the article, in sections 2, 3 and 4, where we 
introduce the parabolic flow on the collision manifold for the 2-dihedral four-body problem, we 
prove the existence or non-existence of parabolic trajectories connecting rest points or singular sets, 
and then analyze the consequences on orbits close to total collision via topological transversality. 
Many tedious computations are postponed to the appendix. 



2 Blow up, regularization and total collision manifold 

We set n := {1, ...,n}. Let V := M. d denote the Euclidean space of dimension d and n > 2 
an integer, and denote with the origin in K d . Given n point mass particles m V P^, k € n 
with positions q^., momenta pk and masses € K + , let q,p £ V n be the vectors (qi, . . . ,q n ), 
(pi, ...p n ). Let a > be a given positive real number. We consider the potential function (the 
opposite of the potential energy) defined by 
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If M denotes the moss matrix namely the diagonal nd x nd matrix given by 



M := diag(mi , . . . , m 1? . . . , m n . . . m n ) 

d times 



then Newton equations 



can be written in Hamiltonian form as 
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We will assume that the center of mass and the linear momentum remains at the origin and we 
define 

n n 

Q = {q=(qi,...,q n ) G V n \ ^771^ = 0}, P = {p = ( Pl , . . . ,p n ) £ F"|^p i = 0}. 



For each pair of indexes i,j E n let Ajj denote the collision set of the i-th and j-th particles 



{q € Q\q t = q 3 }. 



For i ^ j, let A = {J t . Ajj be the collision set in Q. The set of 

collision-free configurations is denoted by Q = Q\A. The differential equations (JlJ then determine 
a vectorfield with singularities on V n x V n , or a real analytic vectorfield without singularities on 
(V n \ A) x V n . The vectorfield given by ([!]) is everywhere tangent to Q x P and so this 2d(n — 1) 
dimensional linear subspace is invariant under the flow. We henceforth restrict our attention to 
the flow on the phase space Q x P. Consequently, H is an integral of the system. This means 
that the level sets S e := H^ 1 (e) fl (Q X P) are also invariant under the flow 0. We observe that 
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E e is a real analytic submanifold of (Q — A) x P having dimension 2d(n — 1) — 1 and the flow 
is not complete, however. In fact certain solutions run off in finite time; this happens exactly in 
correspondence to certain initial conditions leading to a collision between two or more particles 
and the corresponding solution curve meet A in finite time. Since we are interested primarily in 
solution curves leading to total collapse and since the center of mass is fixed at the origin, this 
must occur at the origin of Q. 



2.1 McGehee coordinates 

The aim of this section is to recall McGehee coordinates as given in FP08] in order to fix notations. 
Our basic reference will be |FP08j and references therein. 

The equations can be written in polar coordinates by setting the mass norm in V n defined 
for every q € V n as 

\\q\\ 2 = (Mq,q) 
and suitably rescaling the momentum as follows 

p — \\q\\ . s — — , z — pPp with a — 2/3. 

P 

In these coordinates equations ([I]) can be read as 
(2) 



P = (z,s)P 
s' = M~ 1 z - (z, s)s 

z' = P(z,s)z + — (a), 



where the time has been rescaled by dt — p + @dr (that is, — = p ^ -r)] now the energy can be 

dr dt 

written as 

(3) H - ^p~ a (M~ 1 z, z) - p- a U{s) = p- a (^{M- X z,z) - U{s) 

Let k := d(n — 1) and let us consider the projection (q,p) H> (s, z) from the full phase space 
QxPto the reduced space S k ~ 1 x P (which is the trivial M fc -bundle on the shape sphere S k ~ 1 ) 

Qx P ^ S 1 "- 1 x P. 

We define the parabolic manifold as the projection of all zero-energy orbits (or, equivalently, of 
the zero-energy submanifold of Q x P) in S x P, where S = S k ~ 1 , that is 

= {( s ,z) eSxP: -(M -1 *,*) = U(s)} C S k - X x P. 

The next change of coordinates, due to McGehee |McG74| . is needed for defining the Sundman- 
Lyapunov coordinate v and for the regularization of the parabolic manifold & . Let v, w be defined 
by 

(v =(z,s) 

)w = M -1 z — (z,s)s. 

Then z — vMs + Mw and (w, Ms) = (so that w is also tangent to S), and ([2| can be replaced 
by 

' p' = pv 



(4) 



?/ = \\w\\ 2 + /3v 2 -aU(s) 
s' = w 

<W ' = - |H| 2 s + (p - l)vw + M^WsUis), 
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where V s denotes covariant derivative associated to the Levi-Civita connection induced by the 
Ricmannian metric (-,-)m : — (M-,-), i.e. the component of the gradient tangent to the inertia 
ellipsoid ||<jr|| = 1, i.e. 

f\TT 

V S U = —{s)+aU{s)Ms. 

In fact 

dU 

(M- l V a U{s),v) = ( — (s),v)+aU{s)(Ms,v) = 
3U 

= ( — (s),v)=dU(s){v} = (M- 1 y s U(s),v) M 
oq 

The energy relation becomes 

2p a H =\\w\\ 2 + \\v\\ 2 -2U{s) 
while the parabolic manifold & is then defined by the equation 

v 2 + \\w\\ 2 = 2U(s). 

We observe that differential equation Q gives a real analytic vectorfield with singularities on the 
manifold with boundary := (0,oo) x R x TS; equivalently it defines a real analytic vectorfield 
without singularities on [0, oo) x E x T(S \ A). 

In these coordinates the assumption about the center of mass gives: 

(5) 

n n 

S = {s = . . . , s n ) £ V n \ ™>iSi = 0}, W = {w = (tei, ...,w n )eV n \J2 miWi = 0}. 

i=l i=l 

Let M(h) denote the subset of (p, s,w) defined as follows: 

(6) M(h) := {(p, s, w) G [0, oo) x § x W \ ([§ - ^ hold} . 

Let M+(h) and Mq denote the subsets satisfying p > and p = respectively. All three are 
invariant submanifolds for the vectorfield Note A^o is independent of h. We shall refer to Mq as 
McGehee total collision manifold. We observe also that the parabolic manifold & introduced above 
actually is the projection of McGehee total collision manifold (see |McG74| IDev80| IMoe8T| IMoe83j ) : 
the manifold of (s, w) here is not considered as embedded in the space of (p, s, w) with p = 0. 

Remark 2.1 As already remarked by several authors the effect of these transformations and of 
the time scaling has the effect of gluing a boundary given by p = onto the phase space. Since 
p' = when p — 0, this boundary is invariant under the flow generated by Q. 

By the second equation in Q can be deduced the well-known fact that for < a < 2, v is a 
Lyapunov function on the flow in the parabolic manifold, and therefore the flow is dissipative 
(gradient-like). Moreover, the equilibrium points in Q are the projections of the equilibrium 
points of ([2]) (and the projection is one-to-one in the parabolic manifold), which can be found as 
solutions of 



(7) 

Hence all equilibrium points belong to the parabolic manifold & '. 
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2.2 The linearized flow 

As already observed in section before the vectorficld Q has the equilibrium points given in ([7]); 
moreover as already observed all these points lie on the total collision manifold. We now turn our 
attention to the calculation of the characteristic exponents of the various equilibrium points in 
Mq. Accordingly, let s be a central configuration and let v = \/2U (s). Then using 

(8) 



P = pv 
s' = w 



w 



)\\ 2 s + (p - l)vw + M- 1 V.t7(«) > 



one computes the characteristic exponents for the flow in M(h). To study the eigenvalues of 
the linearized vectorfield at a restpoint 7 = (0, s, 0), we introduce coordinates (Cp>Cs>Ctu) in 
the tangent space to R 2nd+1 . Let g be the Riemannian manifold induced by K 2nd+1 on the 

[2d(n — 1) — l]-dimensional manifold Mo, D and — be respectively the Levi-Civita connection 
and the covariant derivative associated to the Riemannian metric g. The linearization along an 
orbit 70 is represented in local chart by the linear autonomous system 



(9) 



D 

dl 



C(t) = CC(t), Vte 



where ( € 7* [TM (h)j is the vector field along the curve 7 and where the variational matrix A is 
the block matrix represented by 



(l 



(10) 



C := 



\ 



\0 M~ 1 DV a U{s) {P-l)vj 



Devaney has observed that one can guess eigenvector for problem of this type. In fact, given the 
matrix below 



(11) 



B = 







,M- 1 L>V s C/(s) {P-l)v, 



and assuming that £ is a [d(n— 1) — l]-eigenvector of the matrix A := M~ 1 DV s U(s) (the Hessian 
D 2 [/(s))associated to the eigenvalue A, then the [2d(n — 1) — l]-dimensional vector £ = (£, fi£) T 
satisfies the following: 



G8-l)w 



By a straightforward calculation it follows that /i is an eigenvalue of B if 

fi 2 + (1 - P)vpi = A, 

namely 





Ml,2 



(j3 - l)v ± y/(j3 - l) 2 v 2 + 4A 



We observe that, in the Newtonian case corresponding to j3 = 1/2 this formula agrees with that 
of |DV99j . By this we have the following result. 



Proposition 2.2 All equilibrium points are hyperbolic. 
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Proof. For the proof of this result, see |FP08( Proposition 3.8] 



□ 



Proposition 2.3 The dimension of the stable (unstable) manifold of(v,s,0) withv = \/2U(s) > 
is 3 (2) ij s is a rectangle; it is 2 (3) ifsis a tetrahedron. The dimension of the stable (unstable) 
manifold of the point (— v,s, 0) with v — \/2U (s) > is equal to the dimension of the unstable 
(stable) manifold of (v, s, 0). The intersection of the stable (unstable) manifold o/(v, s, 0) with 
the parabolic manifold £P has codimension 0(l)in£Pifv>Q. It has codimension 1 (0) in if 
v < 0. 

Proof. For the proof of this result, see |FP08[ Proposition 3.9] □ 



2.3 McGehee coordinates for the (anisotropic) Kepler problem 

Consider the classical anisotropic Kepler problem, in which S C K 2 is a 1-dimensional manifold. 
Let (0, 2tt) — > S be a (maybe partially defined) local parametrisation, which we will denote by 
9 H> s(6), where s(9) = (m 1 ^ 2 cos 9, m 2 1 ^ 2 sin(9), if mi and m-i are the diagonal entries of M. If 
we introduce the vector s = (—m 1 sin 9, m 2 cos 9) and the scalar w such that w — ws, then 
equations Q turn out to become 

r v' = w 2 + (3 (v 2 - 2U(9)) 

(12) J e' = w 

[w 1 = (p-l)vw + U'{8), 

where we use the primes for denoting the differentiation with respect to r except in the potential; 
moreover the parabolic manifold & defined by the equation 

v 2 + w 2 = 2U{9). 

Now consider the flow in part of the parabolic manifold £P contained in the half-space w > 0: by 
eliminating the term w 2 in the equations of v' and 9' , the projection of the flow on the {9, i>)-plane 
is contained in the region {(9,v) : v 2 < 2U(9)} and is given by the system 



(13) 



which can be written also as 



9' = y / 2U(9) - v 2 

v' = {\-(3){2U{9)-v 2 ) 



dv 



(14) - = (1 - P)y/2U{B) - v*. 



For the projection of the part in w < 0, the first equation of ( 13 ) has to be changed in 9' = 
— yj2U(9) — v 2 . If U(9) = U is constant, that is if the problem is rotationally symmetric, then 
the parabolic manifold is the cylinder of equation v 2 + w 2 = 2U and the flow in 3? (which is 
invariant up to translation in 9) is given by curves leaving the line v = —y/2U, w = of equilibrium 

points at 9 and reaching the equilibrium line v = V2U, w = at 9\ = # H ^, since 

1 p 

(1-(3){9 1 -9 Q )= (2U - v 2 )' 1 ' 2 dv = TT. 

J-V2U 

Thus a "bouncing" trajectory on a collision can be seen as a solution of the regularized problem 

only for 1 — j3 = — implies that a — 1 (this is the reason the Levi-Civita regularization might 

work only in the case a — 1). Nonetheless, the Sundman-McGehee regularization can be applied 
for every a (yielding a less natural regularization). 



G 



2.4 The dihedral 4-body problem 

Let I 3 = C x 1 be endowed with coordinates (z,y), z € C, y € R. For ^ > 1, let Q denote 
the primitive root of unity Q = e 27 ™/ ; ; the dihedral group D[ C 5*0(3) is the group of order 2Z 
generated by the rotations 

Q- ( Z ,V) ^ and k: (z,y) H- (z, -y), 

where z is the complex conjugate of z. The non-trivial elements of Di = (Q , k) are the I— 1 rotations 
around the i-gonal axis Cf i J = 1, — 1 and the Z rotations of angle it around the I digonal axes 
orthogonal to the Z-gonal axis The Newtonian potential for the n-body problem, homogeneous with 
degree —a induces by restriction on the fixed subspace (Ml 6 ') 1 = R 3 a homogeneous potential 
defined for each q £ R 3 by 

as) u{ q )= J2 \i-9ir, 

provided we assume (without loss of generality) all masses mf = 1/1. Now, the potential U in ( [l5| 
can be re-written in terms of coordinates q — (z, y) E C x R as follows. 

On the unit sphere 5* C R 3 (of equation |z| 2 +j/ 2 = 1), parameterized by ((p,9) G (— 7r/2,7r/2) x 
[0, 27r) with y = smip and z = cosipe , the (reduced to the 2-sphere) potential reads 



(16) U (9, ip) = (2 cos ip)~ 



1 



(cos 2 9 + tan 2 ip) a / 2 (sin 2 9 + tan 2 ip) a / 2 



(See [FP08I for the potential in the general case). In spherical coordinates, the symmetry reflections 
of U are (up to conjugacy) 

1. the reflection on the horizontal plane: h v : (9, ip) i-> (9, —ip), 

2. the reflection on the vertical plane hg : (9, <p) (—6, ip) 

3. and the reflection on the plane containing the digonal axis and the point (e OT / 2 ,0), defined 
as h' e : (9, <p) n- (tt/2 - 9, ip). 

Thus we can study U only in the left- upper area of the /^-fundamental domain on S 2 C I 3 , i.e. 
in the geodesic triangle of the shape sphere parameterized by (p,9) € [0,7r/2] x [0,7r/2]. Since in 
the four body case both the reflections h v ,hg give arise the same potential we only distinguish 
two cases. 

1. (Planar case) /i^-reflection (resp. /i^-reflection) : 



2 a [ sin Q 6> cos a 6»J 2 a { sin" ip cos" (p 

2. (Tetrahedral case) (9 = tt/4) 

(18) tfrafoO 



(2cos</?) Q (1 + sin 2 ^)"/ 2 



Remark 2.4 H^e observe that the potential in the planar Newtonian case with 1 = 2 is 1/4 the 
potential of the problem studied by the authors in iSL80j : moreover the potential in the anti-prism 
Newtonian case with 1 = 2 is 1/4 the potential of the problem studied by the authors in JDV991. 
This comes from the fact that we normalized all masses with the value l/l. 

1 We observe that in the general dihedral n-body problem case, the reflections h v and hg induce two different 
potential functions. In fact in this case we have three different class of central configurations: 2i-gon, prism and 
antiprism type. For further details, we refer to |FP08| . 
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V 


K 3 

s a 


9> 

c c 

a a 


rectangular 


> 
< 


3 2 
2 3 


3 1 
1 3 


tethraedral 


> 
< 


2 3 

3 2 


2 2 
2 2 



Table 1: Dimensions of stable and unstable manifolds. 



2.5 Central configurations of the dihedral four body problem with ho- 
mogeneous potential 

As particular case of the results proven in }FP08| . we recall that in this case we have 20 cen- 
tral configurations; in particular 12 of all of these are rectangular and the remaining eight are 
tetrahedral. 

Lemma 2.5 (Planar central configurations) For any a € (0,2) there are exactly 



4 central configurations which are h^-symmetric, and they are on the vertices ( e (2k+lW4 q) 
of the regular rectangle, for k € 3. 



• 8 central configurations which are hg -symmetric (up to conjugacy), and they are precisely on 
the vertices of a prism: (l/v^e^ 7 ™/ 2 , ±l/y/2). 

It is left to compute critical points for 8 = (2k + 1)tt/4, that is, to find fog-symmetric central 
configurations. 

Lemma 2.6 (Tetrahedral central configurations) For any a e (0,2) there are exactly 8 
central configurations which are h' e symmetric (up to conjugacy). They are on the vertices of a 
prism: 

(cos<^e (2fc+1)7rl/4 ,±sin<£). 

Since there are no other central configurations, we can summarize the results in the following 
proposition. 

Proposition 2.7 All central configurations in the dihedral A-body problem are symmetric for one 
of the three types of reflections h v , hg (cfr. Lemma 2.5) or h' g (cfr. Lemma 2.6). 



2.6 The linearized flow 

As already observed in section before the vectorficld Q has the equilibrium points given in ([7| 
and all these points lie on the total collision manifold. Moreover we 

We now turn our attention to the calculation of the characteristic exponents of the various 
equilibrium points in Afo- Accordingly, let s be a central configuration and let v — y/2U(s). 
Then using ^ one computes the characteristic exponents for the flow in M(h). In fact , 

Now we use the above results to describe the set of orbits which begin or end in quadruple 
collision. Orbits which begin at collision are called ejection orbits; orbits which end at triple 
collision are called collision orbits; and orbits which do both are known as ejection- collision orbits. 
It is well-known result that any such orbit must be asymptotic to one of the equilibria associated to 
a central configuration. That is, they lie on the stable and unstable manifolds of these equilibria. 
We denote by E(c.c) and C(c.c) the set of ejection and collision orbits at the central configuration 
c.c As a direct consequences of the dimension of the stable and unstable manifolds, we have the 
following: 
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Proposition 2.8 In any energy level H — h, both C and E consist of union of twenty sub- 
manifolds, where twelve (corresponding to the planar configurations) are bi- dimensional and the 
others are three-dimensional (the ones corresponding to tetrahedral central configurations). All 
ejection orbits emanate from the equilibria v > , whereas all collision orbits are asymptotic to 
the equilibria with v < 0. 



2.7 Regularization of double collisions 

It is not know if a "^-regularization of simultaneous binary collision is possible in general, but in 
this problem, due to the symmetry, simultaneous double collisions can be regularized using the 
type of transformation used by Devaney in |Dev 80 . The aim of this section is to provide the 
regularization of double singularities both in the planar and tetrahedral cases. 



Regularization in the planar case 

We define in the planar and prism case the regularized potential as follows 

Wpi(e) := sin"(20)[/ P1 (0) 

and we consider the new variables 

d sin" (26>) d 



sin" (260 
u := — ; w. 

Then the equation of motions become 
, sm a {29)pv 



do- ' y/Wpi(9) dr' 



(19) 



sin" (20) 



Wm(0)t 



sin 2 " (26») 



+ P[v'-2- 



Wpi(0) 



sin a (26>)03 - l)vu W Pl (8 



y/Wpi(e) 2Wp X {6) 



sin" (26») 



(2 sin" (2(9) - u 2 



2a 



tan(20) 



(sin"(20)-w 2 ), 



where we still use the primes for denoting differentiation with respect to a except in W F ^(9). 
Moreover in these new coordinates the energy relation becomes 



(20) 



.2 , ,2 o„a m sin 2Q (20) 



u z + {v z -2p a H)- 



Wpy{6) 



2sin"(260, 



and by taking into account of the energy relation, the equation of motions could be written as 

follows 

(21) 

, sin a (20)yOT 



yWp7(0) 

V' = y/W P1 (0) 



sin" (20) (/3- l)vu W } 



Piy " (2sin"(20)-u 2 ) 



2a 



(2p a H- w 2 )sin"(20) 



y/W Pl (0) 2Wpx{0y v ' ' tan(20) 

(We recall that the Newtonian case corresponds to a = 1). Explicitly 

/sin" (20) sin". 



Wpi{6) = sin" cos" + sin" + cos" 0, W Pl {6) = a 



tan(20) tan0 



cos tan ( 
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Remark 2.9 There is a strict analogy between these equation of motions and the differential 
system given in \Dev8(k Eqn.(1.12)] and in WV991 Eqn. (2.3)]. We also observe that the above 
regularization holds equally well for the parabolic manifold. In fact, the variable p is not essential 
neither to regularize the total collision manifold nor to investigate the flow on the total collision 
manifold. 

By taking into account the energy relation and by putting p = 0, it follows 
(22) f+ji^^W, 



Regularization in the tetrahedral case 

We define in the tetrahedral case the regularized potential as follows 

W To t{p) := (2 cos p) a U Tct {p), 
and we consider the new variables 



(2cos^) Q 
u := — : w, 

\/W T ct{p) 



d _ (2cos(p) a d 



Then the equation of motions become 



(23) 



(2 cos p>) a pv 

yJW Tet {p) 

(2 cos p) a 



W Tet {p)u 2 

(2 cos p) 2a 



+ (3 v 



, W T et(y) 
'(2cOS(,9) Q 



p = u 



(2cosp) a (B — l)vu W^ nt (p r „,„ .„ r ,„ .„ 9l 

Tetv ^ [2(2 cos ip) a -it 2 ] +atantp[(2cos(p) Q -it 2 ], 



y/W T ct(p) 2W^ctM 

where we still use the primes for denoting differentiation with respect to a except in W^ ot (<p). 
Moreover in these new coordinates the energy relation becomes 



(24) 



i 2 + (v 2 2p»H) = 2(2cob y)", 



and by taking into account of the energy relation, the equation of motions could be written as 
follows 



(25) 



, (2 cos (p) a pv 
p = — , 

\/W Te t{<P) 

v' = y/W Te t(p) 



W Tei (p>) 
tp' = u 



+ 



[2(2 cos tp) a - u 2 } + 



\/WrM 2W Tc t(p) 

(2p a H - v 2 ))(2cosp) a 



Explicitly 



W Tet (p) = 1 + 



a tan p>(2 cos p) a 
2 i+P cos « v 



W Te t(<p) 



(1 + sin 2 pf ' 



Wn* = - 



a sin p cos" 1 p 



(1 + sin 2 p)P+ l 
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Remark 2.10 There is a strict analogy between these equation of motions and the differential 
system given in \Dev8(k Eqn.(1.12)] and in WV991 Eqn. (2.3)]. We also observe that the above 
regularization holds equally well for the parabolic manifold. In fact, the variable p is not essential 
neither to regularize the total collision manifold nor to investigate the flow on the total collision 
manifold. 



By taking into account the energy relation (24 1, the regularized total collision manifold reduces to 
(26) u2+ ^ 2 (2cos^ = og 

We shall discuss the properties of this flow in the next section. 



3 Colliding and non-colliding parabolic connections 

The object of this section is to describe in details the flow generated by the differential system 



(191 and (25) on the total collision manifold in the planar and tetrahedral case, respectively. A 
direct consequence of this study is the proof of the existence of connecting orbits on the invariant 
subset of the 3? fixed by the the reflections hg,h,p,h' s . The idea to perform our analysis is based 
upon a study of the intersection between the stable and unstable manifold of the equilibria on the 
parabolic manifold 3 s . 

Due to the gradient-like character of the (regularised) flow on 3 21 , each orbit in C 3 s either 
tends to a rest-point or v — > +oo as t — > +oo. Moreover by taking into account the inequality 
v 2 < 2U(s), it is clear that if v —> +oo then U(s) —> +oo and so two or more particles have to 
collide. Thus if s is parameterized in cartesian coordinates by (si, S2, S3), this means that one of 
the distance sj + s 2 — > for 1 < i < j < 3, that is the particles tend to have a binary configuration. 
Therefore, by following |Vid991 Section 6], we can define the stable binary escape sets as follows: 

(27) := {(«,«, to) G &>: v — > +00 and Sj — > ±1 for t — > +00} , 
for j G 3. In an analogous way we define 

(28) ^" ,=fc := {(v,s,w) G 9>: v -> -00 and Sj ±1 for t -t +00} . 

The binary sets together with the 40 critical points are the only possible asymptotic alpha and 
omega limits for an orbit on 3P . 

Notation 3.1 If X and Y are two equilibrium points or even two binary escape sets on & we 
shall write X Y in order to indicate the existence of an orbit in £P with asymptotic behavior 
X for t — > —00 and Y as t — > +00 i.e. if there exists an orbit on the total collision manifold 
which lies in W U (X) n W S (Y). Since at each central configuration corresponds two value of the 
function v, we introduce the following convention. We denote with superscript + (resp. —), central 
configurations corresponding to the positive (resp. negative) value of the coordinate v. 

Remark 3.2 In order to simplify the study of the connecting orbits between critical points of the 
function v some remarks are in order. 

1. At first we observe that the presence of the symmetry of the vector field Q, namely 

(v, s, to) n- (— v, s, —w) 

implies that for each connection X + Y + there is a symmetric one Y~ X~ which will be 
term dual connection. Moreover this transformation sends the stable and unstable manifold 
of the point X + respectively on the unstable and stable manifold of the point X~ . 

In conclusion the transformation (s,w) 1— > (s, — to), together with a time reversal t t— > —t, 
takes orbits into orbits. This transformation carries 3$ s ^ to 38^ and viceversa. 
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The second observation is that the increasing character of the function v imposes some re- 
striction about the existence of connecting orbits between central configurations. In fact, since 
v is non- decreasing and strictly increasing away from the central configurations, connecting 
orbits between X and Y cannot occur if v(Y) < v(X). Denoting by PI, Tet respectively the 
planar and tetrahedral type central configurations, as already observed, we have 

u(Pl) > u(Tet) > 0. 

In fact, in the Newtonian case (a = 1) the inequality above, readily follows since 



v+(V\) := y/l + 2V2 > := «+(Tet). 

However in the homogeneous case, it is enough to observe that the function I : [0, 2) — > M 
defined by 

1(a) := v + (Pl) - v + (Tct) = \J^-°- + 2 2 -"/ 2 - \/2 1 - a Z a Q- a / 2 + 4i-«3 Q / 2 2 Q / 2 , 

is a non-negative increasing function such that 1(0) = 0; thus it is strictly positive for a^O. 

In order to study the possible connections in the dihedral four body problem for any homogenous 
potential of degree —a, for a € (0,2) we shall analyze two cases that it contains as subproblems. 

• The Planar problem. This particular case appears when we consider the invariant set {<ys' = 
ip = 0} (or {8 = 0, 9' = 0}) which corresponds to the subset 3? vx of the parabolic manifold 
fixed by the symmetry h v . In this case the parabolic manifold is two-dimensional and in 



the Newtonian case a = 1 it was studied by }SL82j : in Figure 3.1 it is given a qualitative 



description of the flow on the covering of the regularized component. 

The Tetrahedral problem. This particular case appears when we consider the invariant set 
{8 = ir/4,6' = 0}. Also in this case the parabolic manifold is also two-dimensional. This 
problem in the Newtonian case was studied by Delgado & Vidal in [DV99 , and the flow on 
the covering of the regularized component is qualitative depicted in Figure [2] 

The main goals of this section is to completely study the flow in the above two cases. In Theorem 
|3.11| are described all connecting orbits corresponding to the planar section and in Theorem |3.13| 
and Theorem |3.16| are described all connecting orbits corresponding to the tetrahedral section. 
Before proceeding further with the description of the flow on the parabolic manifold some com- 



ments are in order. First of all in figure 3.1 it is represented the flow obtained by integrating the 
corresponding 1-dimensional Cauchy problem for different values of the initial conditions. This 
flow is represented on the covering over each piece of 1-sphere between the binary collision and 
the binary collision. Figure |3.1| actually represents the flow on the regularized component of the 
parabolic manifold over the 1-sphere corresponding to the planar section. Figure [2] represents 
the flow on the covering of a regularized component of the parabolic manifold for the tetrahedral 
section. In this case the situation is more complicated due to the presence of two different types of 
central configurations. By changing coordinates, locally in the neighborhood of binary collisions, 
the regularized component of the parabolic manifold in the planar section resemble that studied 
by McGehee in |McG74] for the collinear three body problem, while the regularized component 
of the parabolic manifold in the tetrahedral resemble that studied by Devaney in [Dev80j in the 
isosceles three body problem. 
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3.1 Flow on the parabolic equation: the planar case 

The projection of the equation of motions on the parabolic manifold are 

v 2 sin" (20)" 



(29) 



v' = (1 - PWWpi{6) 2- 

6' = u 
, _ sin a (29){f3-l)vu 



W P1 ( 



(2 sin" (26>) - u 2 ) 



2a 



y/Wp~$) 2Wpi{6) 

where the regularized parabolic manifold 2? is given by 



tan(20) 



1 - 



2 sin Q (26») 
Wpi(5) 



(30) 



2 2 sin 2 " (2(9) 



= 2sin a (26»), 



Now consider the flow in part of the parabolic manifold 3? contained in the half-space u > 0: the 
projection of the flow on the (9, «)-plane is contained in the region {(6,v) : v 2 < 2Up\(8)} and is 
given by the system 

(31) 

which can be written also as 
(32) 



9' = ^2U Pl (9) - v 2 

v' = (l-f3)(2U P1 (&)-v 2 ), 



dv 



= (l-P)y/2U m (d)-V*. 



For the projection of the part in u < 0, the first equation of (131 has to be changed in 9' = 
-^/2U{9)-v 2 . 

Notation 3.3 In order to give the complete description of the flow on the parabolic manifold 
including orbits that run out along the arms of binary collisions, we use the following convention. 
Let W s (jr/2) (W u (tt/2)) denote the set of points on whose forward (backward) orbit run out 
along the upper (lower) arm at 9 = n/2. Similar definitions for W s (— tt/2) (W u (—tt/2) and W s (0) 
(W u (0)). 

Before stating the main result of this section, we remark that the projected differential system 
(36) has the following symmetries: 



(33) 



and the composition of both. 



(\ : {9, v, u, a) i — > (9, —v, —it, — a) 
$2 : (9, v, u, a) i — > (tt/2 — 9, v, —u, a) 



Proposition 3.4 For any a € (0,1], the global flow in the parabolic manifold for the dihedral 
planar four body problem is described by the following relationship among the stable and unstable 
manifolds on 2P: 

1. W u {c_) C W s (0) U W s (-k/2). 

2. W u {eir/2) n W s (cm/2) is an open set for e, a € {0, 1}, 

where we denoted by c the planar central configuration corresponding to the negative value of v. 
Proof. The proof of this proposition is deferred to the Appendix [Aj □ 
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Figure 1: The flow on the covering of a regularised component of the parabolic manifold for the 
planar section. 



Remark 3.5 Proposition 3.4 implies that there are no saddle connections on , thus preventing 
regularizability of the total collision. 

In our case item 1 in Proposition \3.4\ says that orbits passing close to the planar homothetic 
orbit escape from a neighborhood of total collapse with successive binary collisions; item 2 says 
that any combination of orbits running from one arm of binary collision to another can occur, and 
there is an open set of initial conditions in 3? whose orbit have this property. In particular no 
saddle connections on £P , thus preventing regularizability of the total collision. 

Proposition 3.6 For any a G (1,2), the global flow in the parabolic manifold for the dihedral 
planar four body problem is described by the following relationship among the stable and unstable 
manifolds on . There exist ao*,a, <E (1,2) with 1 < ao* < a* < 2, such that 

1. W u (c-)cW s (0)UW s (ir/2),foranyae (1, 2) \ {a *, a*}. 

2. For a = ao* the right-hand branch of W u (c_) coincides with the right-hand one of W s (c+). 

3. For a — a* the right-hand branch of W u {c_) coincides with the left-hand one of W s {c + ). 

4. W u {en/2) n W s {o-k/2) is an open set for e, a € {0, 1}, 

where as before we denoted by c^(resp.c + ) the planar central configuration corresponding to the 
negative (resp. positive) value of v. 

Proof. The proof of this proposition is deferred to the Appendix [Aj □ 



3.2 Flow on the parabolic equation: the tetrahedral case 



The projection of the equation of motions on the parabolic manifold and in the tetrahedral case 

are given by 

(34) 

u 2 (2cos<p) Q 



v 1 = {1-[3)^WtM 



W Tet {ip) 



ip = u 

, (2cos^H/3-l) TO WiaJtp 



[2(2 cos tp) a - u 2 ] - a tan <p(2 cos ip) a 



1 - 



v 2 (2 cos tp) a 
W Tc t((p) 



14 



and the regularized parabolic manifold takes the form 

2 , 2 (2C0S(^) 2Q 



(35) 



W T etM 



2(2 cos ip) a . 



Now consider the flow in part of the parabolic manifold & contained in the half-space u > 0: 



the projection of the flow on the 
is given by the system 

(36) 

which can be written also as 
(37) 



, v)-plane is contained in the region {{6, v) : v < 2Up\(9)} and 



v' = (l-P)(2U Tei (<p)-v 2 ), 



dv 
dtp 



= (l-/3)\^p 



For the projection of the part in u < 0, the first equation of ( 13 ) has to be changed in ip 1 
-^2U Tet {v)-v 2 - 



The differential system (361 has the following symmetries: 



(38) 



S\ ■ (p, v, u, a) 
y 2 ■ (<P, v, u, a) 



((p,-v,-u,-a) 
(-(f,v,-u,a) 



and the composition of both. The first symmetry comes from the reversibility while the second 
comes from the reflection with respect to the horizontal plane. 

Proposition 3.7 The global flow in the parabolic manifold for the dihedral planar four body prob- 
lem is described by the following relationship among the stable and unstable manifolds on 

1. W u {c_) C W s (-tt/2) U W s (tt/2). 

2. W u {eir/2) n W s (cm/2) is an open set for e, a £ {0, 1}. 

3. W s (r_) C W u (c_) U W u {tt/2) 

where we denoted by c_ the planar central configuration corresponding to the negative value of v 
and r_ is the tetrahedral central configuration. 



Proof. The proof of this proposition is deferred to the Appendix [B] 



□ 



Remark 3.8 In our case item 1 in Proposition \ 3. 7| says that orbits passing close to the planar 
homothetic orbit escape from a neighborhood of total collapse with successive binary collisions; 
item 2 says that any combination of orbits running from one arm of binary collision to another 
can occur, and there is an open set of initial conditions in & whose orbit have this property. The 
last item says that the stable branch of r_ die in c_ and in the the lower arm of a double collision 
ate = it/2. 

Proposition 3.9 For any a 6 (1,2), the global flow in the parabolic manifold for the dihedral 
planar four body problem is described by the following relationship among the stable and unstable 
manifolds on & . There exist ao*,a* € (1, 2) with 1 < «o* < a* < 2, such that 

1. W u (c_) C W s (-it/2)UW s (it/2), for any a e (1, 2) \ {a *, a*}. 

2. For a = ao* the right-hand branch ofW u {r_) coincides with the right-hand one ofW s {l+). 

3. For a = a* the right-hand branch of W w (r^) coincides with the left-hand one ofW s (r + ). 
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Figure 2: The flow on the covering of a regularised component of the parabolic manifold for the 
tetrahedral section. 

4. W u (en/2)nW s ((nr/2) is an open set for e, a G {-1,1}, 

where as before we denoted by C-(resp.c+) the planar central configuration corresponding to the 
negative (resp. positive) value ofv, r_(resp.r + ) the tetrahedral central configuration corresponding 
to the negative (resp. positive) value ofv for ip > and finally l-(resp.l + ) the tetrahedral central 
configuration corresponding to the negative (resp. positive) value of v for ip < 0. 

Proof. The proof of this proposition is deferred to the Appendix [5] □ 

Notation 3.10 We denote the central configurations in the first octant of the shape sphere cor- 
responding to the fundamental domain by: 

1. (tetrahedral c.c): e/-i := (A;7r/4, # c ), e/-2 '■= {kit/A, — 6 C ) 

2. (planar c.c.): p lk := (kn/^O), p 2k ■= (n/2,n/4), p 3fc := (0,tt/4) 
where 9 C := arctan(l/v / 2) and k G 3. 

Theorem 3.11 There exist ao,*, a* G (1, 2) with ao* < a* such that on the parabolic manifold & 
the following connections and their dual ones occur: 

1. for a G (0,a) \ {a *} 

(PI) p+-^< + ; 
(P2) p+^^ s 2 ' + ; 

(PS) ^< + ^ P ti; 
(P4) ^' + ^p!i- 

(P5) There is no saddle connections between two planar central configurations PY^ . 
Moreover in the case (P4) and (P5) the last binary collision occur for a positive value ofv. 

2. for a = ao* there are saddle connections p^ ^ pj^ and p^ ~> p\ 2 passing through one 
binary collision. 

3. For a = a* there are saddle connections p^ pf x , p^ ~> pf 2 and finally p^ ~+ pf 3 
passing through two binary collisions. 
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4- For a G (a*, 2) 

(P3) az'+^ptti 
(P4) £%< + ^ P i v 

However in this last two cases the last binary collision occur for a negative value of v. 

Remark 3.12 We observe that analogous connections and saddle connections occur between the 
remaining central configurations and (or) central configurations/binary escape sets. 

Proof. The proof of this result immediately follows by the discussion at the beginning of this 
section as well as Proposition |3.4| and Proposition |3.6| 



Theorem 3.13 There exist 70,*, 7* € (1, 2) with 70* < 7* such that on the parabolic manifold £P 
the following connections and their dual ones occur: 

1. for^e {0,~f*)\{lo*} 

(Tl) ejj ~> Pxx, ~> P11 (non- colliding parabolic orbit); 

(T2) e+^@ 3 s ' + ,e+^^-; 

(T3) ^'-^e+,^'-^e+; 

(T4) %T - e+ . 

(T5) There is no saddle connections between two central configurations. 

2. for 7 = 70* there are saddle connections ^12? ^11 ^21 and finally ^31 passing 
through two binary collision. 

3. For 7 = 7* t/iere are saddle connections ~> and e7 x ~> passing through one 
binary collision. 

4. For a £ (7*, 2) 

(T6) ^-^e+ l! ^-^e+; 
(T7j ~» e+. 

However in this last two cases the last binary collision occur for a negative value of v. 

Proof. The proof of this result immediately follows by the discussion at the beginning of this 
section as well as Proposition |3.7| and Proposition |3.9| □ 



Remark 3.14 Analogous connections holds between the remaining six tetrahedral type central con- 
figurations and the remaining three planar central configurations located on the equator if = of 
the shape sphere. 



3.3 Non-colliding parabolic connections 

The scope of this section is to prove the existence of non-colliding parabolic connections. This will 
be achieved by using topological arguments similar to those used in [Moe83J. Before, we recall 
some basic definitions about topological transversality . 

Given two submanifolds _P, Q of the manifold AT, we say that their are transverse at p if 
pePnQ and if T p P + T p Q = TN. Of course by the rank theorem it follows that the sum of 
the dimension of P and of Q can be greater than or equal to the dimension of the manifold TV. 
Moreover we say that P meets Q transversely if either P D Q = or else P is transverse to Q 
at each p £ P n Q. The important fact is that the transverse intersection cannot be removed by 
a c i 1 - small change; thus in other words, two submanifolds which meet transversely do so stably. 
This statement can be easily proven since the transversality is an open condition. 
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Let P be a p-dimcnsional submanifold of the n-dimcnsional manifold N and let us denote by 
D k the fc-dimensional disk. Given a continuous map of pairs a: (D n ~ p , D n ~ p \{0}) —> (AT, N\P) 
such that c(0) = p G P, we will say that it is a (n — p)- complementary cell at p to P if in a small 
enough neighborhood J7 of the point p we have H n _ p (U, U\(U (1 P)) = Z and er is a generator 
of the (n — p)-dimensional relative homology group. Now we are in position to state the following 
useful definition. 

Definition 3.15 LetP,Q be a p and q- dimensional submanifold respectively of 'then- dimensional 
manifold N . We say that they are topologically transverse at the point p G P n Q and we write 
P ft) Q if there exists a (n — p)- complementary cell at p to P in (Q, Q\Q fl P) and a (n — q)- 
complementary cell at p to Q in (P, P\Q (1 P). 

Theorem 3.16 The following non-colliding connections and their dual ones occur on the parabolic 
manifold & : 

( Ql) e£ Pn, e+ p+, e+ ~> p+; 
fQ2j e+ ^ e+ ~> e+ <v ^ 3 S <+ 

Furthermore for each of these connections the intersection between the stable and unstable manifold 
is topologically transverse. 

Remark 3.17 We observe that analogous connections occur between the remaining central con- 
figurations. Item (Ql) was already proved in Theorem 3.13 Item (Tl) in complete different ways. 

Proof. We shall denote by the symbol &\ the level set of v, namely 

0> x : = {( v ,s,w) G &\ v> A}. 

Clearly since v is gradient-like, it follows that for all A G K each of these sets are invariant for the 
flow induced by the function v. 

Let p: & — > & (defined by (v, s, w) i-> s) be the canonical projection onto the fundamental 
domain of the action of on the shape sphere and we let S\ := p(3P\). Taking into account the 
inequality v 2 < 2U(s) it follows that 

$>\ = {se&: U(s) > A 2 /2}. 

Denoting by c := ^(e^) the level set of the function v at the tetrahedral type central configuration 
then for all positive e the set given by: 

® c+e := {seS>: U(s)>(c + e) 2 /2} 

is just the fundamental domain with a small disks about en removed. The proof of this result is 
completed by using the following lemma. 

Lemma 3.18 Let a: S 1 — > ^ c + e be an analytic closed curve such that p o a: S 1 — > 2> c +e * s n °t 
contractible in @ c +e- Then a intersects and contains one-cells complementary to W s (pf 1 )D 

for each i G 3. 

Proof. The proof of this result follows, up to trivial adaptation, the proof of by [MOC83, Lemma 
4.5 1. □ 



Remark 3.19 Another way to proceed in order to prove the existence of the heteroclinic connec- 
tions predicted in item (Ql) (up to transversality) is to argue as in the proof of \Vid9!^. Theorem 
3]. 
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Figure 3: All colliding and non-colliding connecting orbits on the parabolic manifold with their 
dual ones between the central configurations contained in the fundamental domain and the alpha 
and omega limit sets. 

4 Behavior of orbits close to quadruple collision 

In this section we show how to use the connecting orbits on the in order to establish the motion 
of the dihedral four-body problem. We briefly recall the classical Palis inclination's lemma. 

Lemma 4.1 Let p be a hyperbolic rest-point of a flow on a manifold N. If a positively invariant 
manifold P has a non-empty topologically transverse intersection with W s {p), then W u (p) C P. 

The following result stated by Moeckel in the case of the three body problem, (cfr. [Moc83, 
Proposition] and references therein) has an analogous in the situation we are dealing with. 

Proposition 4.2 Let r,h be given constants, with r > 0. Then there is a constant c(r,h) such 
that any orbit of energy h beginning in the region where p < r and v > c(r, h) is of the following 
type: two opposite sides of the rectangle formed by the four bodies are always the shortest and 
remains bounded while p tends monotonically to infinity. 



Proof. The proof of this result follows by the planar and tetrahedral energy relations ( 20 1 and 



(26) respectively and by the fact that at each instant the bodies form an orbit of the dihedral 
group Z?2- In fact, for large values of v the conservation of the energy forces a small value of the 
multiplicative factor of sm 2a (29)/Wpi(9) in the planar case and of (2 cos <£) 2Q /Wret (</?)• In each 
of these case this condition correspond to the fact that the configuration in which two particles 
are close together. Now, by taking into account the symmetry of the problem, also the remaining 
two are to be closed, and this conclude the proof. □ 



Proposition 4.3 For any j G 3, the following inclusions between invariant manifolds hold, 
(i) W'fcjJcW^j; 
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(ii) W u (p+ 1 ) C ^«(e+). 



Proof. We observe that, by symmetry, it is suffices to show that W u (p~j 1 ) C W u (e'^ 1 ), and this 

will follow from W u {p+ 1 )D & C W u (e+ 1 ) n & , since the subset of the full phase space in McGehee 
coordinates corresponding to {r > 0} is open. By arguing as in the proof of |Moe83[ Proposition 
5.1], the result follows by using Lemma 4.1 and Theorem 3.16 □ 



Corollary 4.4 Arbitrarily close to every planar p^ type collision orbit there are tetrahedral e xl 
type collision orbits. Moreover, arbitrarily close to every tetrahedral type ejection orbit there 
are planar p% type ejection orbits. 

To investigate the behavior of orbits which do not actually collide we consider the connections 
between an asymptotic set for a negative value of v and one for a positive value of v. 

Proposition 4.5 Suppose the following connection occur £ffi' + for i,j 6 3. Then there 

is an open set of orbits with the following behavior: r(t) — > oo as t — > ±oo, but as < — > — oo two 
opposite edges of the rectangle remains bounded (and the remaining ones become unbounded) while 
for t — > +oo the behavior of the opposite edges interchange. 

Proof. Any orbit in M + {h) passing close to an orbit in M$ connecting 3§™ ,+ to enters the 

domain of applicability of Proposition |4.2| in forward and backward time. The conclusion easily 
follows. □ 



Proposition 4.6 Suppose t ^9^' where t is any rectangular or central configuration. Then 
arbitrarily close to every t~ collision orbit in {r > 0} 7 there are orbits with the following behav- 
ior: after reaching a minimum, r(t) tends monotonically to infinity while in a rectangular type 
configuration in which two edges become bounded and the others longer and longer. 

Proof. Let tyf be any neighborhood of a point p on a t~ collision orbit and let us consider the 
forward orbits of the points in a i/ ', namely % ■ [0, oo), the union of the forward orbits of the points 
on % . Sinc e °M ■ [0, oo) is topologically transverse to W s (t~) and positively invariant, by using 
Lemma 



4.1 



we conclude that W u {t~) C tf£ ■ [0, oo) and therefore there are points in ^ which pass 
arbitrarily close to a t connecting orbit. But to any orbit passing close enough to such 

an orbit it is possible to apply Proposition |4.2| and the conclusion easily follows. □ 

We turn now to a consideration of the more delicate phenomena associated to a topologically 
transverse connections between a central configuration corresponding to a negative value of v and 
one corresponding to a positive value. 

Proposition 4.7 Suppose that the topologically transverse connection t u + on the parabolic 
manifold occurs. Let p € W s (t~) and q G W u (u + ) and let ^ p and ^ q be any neighborhoods 
of these points in M + (h). Then there is an orbit segments in M + (h) which begins in °tt v , passes 
close to a triple collision, and ends in °l/q. 

Proof. Since the dimension of the stable and unstable manifolds are as in the case of the spatial 
three body problem, the proof of this result is a trivial adaptation of [Moe83| Proposition 5.5]. □ 



Remark 4.8 Some remarks are in order. First of all we observe that all these results reveal in a 
certain sense a chaotic behavior of the total collision orbits in the sense that all these orbits are 
extremely sensitive to small changes in initial conditions. For instance in the first result of this 
section we proved that after passing close to a planar type configuration, the system emerges with 
arbitrarily large kinetic energy. In particular in every neighborhood of every planar type collision 
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(a) Four bodies with dihedral sym- (b) Four bodies with dihedral sym- 
metry and anti-symmetric con- metry and semi-coreography con- 
straint, straint. 



(c) Twelve bodies with dihedral 
symmetry and 3-coreography con- 
straint. 



orbit there are tetrahedral collision orbits as well as orbits which avoid collision and emerge from 
a neighborhood of the singularity in any of the irregular rectangle configuration. 

Moreover, when a total collapse is reached (asymptotical to a planar central configuration), 
the solution oscillates very rapidly passing through an arbitrarily number of instantaneous planar 
configurations (this is due to the spiralling character of the sink) until a double collision occurs, 
followed by a passage through a planar central configuration and further escape from a neighborhood 
of total collapse with multiple binary collisions. An approach to total collision asymptotically to a 
tetrahedral type central configuration is very unstable in the sense that after a close encounter is 
a double collision followed by a planar configuration and then an escape from total collision with 
multiple binary collisions. 



Remark 4.9 Variational methods under symmetry constraints have been popular in the recent 
years, starting from Chenciner and Montgomery "figure eight" orbit [Che02, CM001. Symmetric 
periodic minimisers have been found by many authors (to name a few, \A CZ94\ I CheOS\ \ Che 03a. 
\Che0 3b. SX05i). Periodic orbits dihedrally symmetric can be shown to exist as well: from JFer06, 
\Fer07\ [FTO^j it follows that ifx(t) is a minimiser of the Lagrangean action functional, constrained 
to the Sobolev space of G-equivariant periodic trajectories (where G is a cyclic extension of Di), 
then x(t) is collision-less. In the Figures \4(a)\ \4(b)\ an d 4( c ) ' some of the resulting orbits are 
shown, which are symmetric with respect to some cyclic extensions of the group Di . 



A Some pointwise estimates for the planar case 



A.l Proof of Proposition 3.4 



The proof of this result will be divided into two main steps. The first correspond to the Newtonian 
case a = 1 and the second to the case a € (0, 1). 

Moreover we refer to the right (left) unstable branch of W u (c_) as that having u > (u < 0) 
in a small neighborhood of c in 3?. 

Using uniqueness of solutions, the symmetries and the gradient-like character of the flow with 
respect to v, it is sufficient to prove the following fact. 

(1) The right branch of W u (c-) intersects 9 = ir/2 (9 = —tt/2, resp.) with v < and 
then intersects the section v = with an angle 9 € (0, 7r/4). 
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Proof in the Newtonian case. 

The proof of this claim will be divided into some parts: 
(la) Upper bound for v(tt/2) in the case u > and 9 6 (tt/4, ir/2). 



dv 



i 



i 



i 



cos 9 sin 9 



< 



< 



1 1 



a 2 



where 



Thus 



a? := max 



v(n/2) < u(7r/4) + 



v(ir/4) 



2 V cos 

f{6) = l + V2, f(9) = 1 + -1 
*/ 2 l / i 



sin I 



1 



tt/4 2 V COS 
tt/4 



a 2 d# =$ :=Tr /2- 



1 



sin$ 



a 2 d$ < 



< «(tt/4) + - 

z Jo 



tt/4 



ai 



a 2 d$ 



where cti := max[o iJr /4] — 



-\/2 since the function <?($) 



4> 



sin$ 4_ sin$ 

[0, 7r/4]. Since v (tt/4) = — v/l + 2\/2 and by using the quadrature 



is increasing on 



*/ — + bdx = y/x(a + bx) H -= log[2\/&\/ (x(a + 6a;) + 26x + a], 

2V6 



it follows that 



v(n/2) < v x := 



1 



tt/4 



'■Ky/2 
Ax 



1 + V2 - V 1 + 2^2 w -0.8014 < 0, 



therefore v intersects = n/2 at a negative value. 

(lb) Lower bound for v(tt/2) in the case u > and # € (7r/4,7r/2). 

Next we calculate a lower estimate for v(ir/2). For this we first observe that 



dv _ 1 / 1 1 
d9 ~ 2\ + shifl + cos" 



z; 2 > 



1 1 



2 V cos( 



+ a 3 - v 2 , 



where 0,3 := f(d) — 2. Now we consider the interval J := [3/8-7T, 7r/2] C 

[7r/4, 7r/2]^]and by a direct calculation: 



|>i^L +2-^/4) ^(*/2)>«(3/*r) + ± 



r/2 



3/87T 



COS( 



+ 2 -w 2 (tt/4) d6»> 



> u(7r/4) + 



1 



4> 



a 4 d<i> 



1.4164 := v 2 < 0, 



for $ as above and 04 := 2 — v (tt/4) = 1 - 2V2; moreover in the last inequality we used 
the fact that sin$ < From (la) &(lb) it follows that v 2 < v{9) < v 1} for 9 e [7r/4,7r/2] 
and u > 0. 



2 We observe that the value 8 := 3/87T has been chosen suitable in order the estimates to work and in such a way 
the expression under square root is nonnegative. 
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(lc) The continuation of the right branch through double collision at 6 = n/2 for u < 0. We 
affirm that v remains negative in the interval J. 

In fact, for u < 

dv i ri 

~dB - 2V <5osl + ° 5 

for a 5 := f (3/8tt) = 1 + 2 1 w 2.0824 and/ as above. Let a 6 := ff(7r/8) = 

V2 + \/2 

3/4 = w 1.0260. Integrating we get: 

V2 + V2 



K3/8.)<^/2) + ^^^ + a 5 ^< 

< v(n/2) + i ^ +a 6 rf$ < Wl + 1 ^ +a 5 d$ 

«Dil 0.7111 « -0.0903 < 

Therefore 

«2 < u(tt/2) < v(6) < v(3/8n) < 0, for 6 £ [3/8tt, tt/2], and u < 0. 
Therefore 

^2 < u(tt/2) < u(0) < u(3/8tt) < 0, for 9 e [3/8tt, tt/2], and u < 0. 

(Id) The continuation of the right branch through double collision at 9 = n/2 for u < 0. We 
affirm that v remains negative at the point 9 — it/A. 

In fact, for u < 

dv i ri 

- 2V <5o7# +ct5 

for a := f(n/A) = 1 + 2\/2 and/ as above. Let b := g{n/A) = n/A\f2. Integrating we get: 

v(tt/4) < v(n/2) + - J —+ad9< 

2 J n / 4 V cos 9 

(40) i r /4 /6 , i r /4 /& 



'3/8 

( 39 ) , , , 1 r /8 , 1 



" v{n/2) + 2 J \l$ +aM - Vl + 2j \$ +aM 
w wi + 0.7445 = -0.5630 < 

(le) From the above inequality we get the estimates v\ > v(9) 2 , for all 9 £ [37r/8, tt/2] and u < 0. 
Using this estimate we can write: 

dv 1 [ 1 1 l/i 

-M = 2y 1 + ^9 + ^T9- V ^2V^ +a7 - V * 



where a*? := min ee [ 3 / 87r 7r /2] /(#) = J{' K I2) = 2. By a direct integration we have: 

i r /2 n 

w(3/8tt) > v(tt/2) + - x + a 7 -v 2 2 d9> 

1 J3/8n V cos 9 

>v(n/2) + -J^ x J- + 2-v 2 d^ 



>v 2 + - a/- + 2-v| rf$ w -0.7900. 
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(If) To finish the estimate for the continuation of the right branch through double collision, we 
now consider £ [0, 37r/8] and u < 0. We observe that in = ir/4, U(9) has a relative 
minimum and therefore for £ [0,3tt/8] we have U{9) > U(w/4). Now let v := ^2[/(tt/4); 
thus 



dv 

~d9 ~ 2 

and by a direct integration, we have 



arcsm i 



(-) > / d9 = 3tt/16. 

Vu/ „(o) 2J 



Therefore 



■u(O) > 5sin(37r/16 + arcsin(u(37r/8)/5)) > 5 sin (3tt/16 + arcsin(-0.7900/£)) w 0.3376. 

We observe that this estimate could be refined by considering the interval [0,7r/4] instead of 
[0,37r/8]; however to do so we need an estimate from below for v{tt/4). 

By symmetry, the same arguments holds also for the left branch of W u (c_ ) . Now the thesis easily 
follows by the fact that the flow is of gradient-type with restpoints c_ and c+. □ 



Proof in the case a £ (0, 1) 

Also in this case the proof of this claim will be divided into some parts: 
(la) Upper bound for v(ir/2) in the case u > and 9 £ (7r/4, 7r/2). 



dv 1 L 
ci0 2\ 



cos" sin" 9 



- v 2 < 



1 



2V2+0 



cos" sin" 9 



< 



< 



1 



2V2+/3 V COS" 



+ b 2 



where 



Thus 



b 2 := max = 1 + 2 P , h{9) := I + 

ee[7r/4,7r/2] sin" 



2V2 

= «(*/4) + ^ / 
< w(tt/4) 



/2+/S 



1 



COS 



— + b 2 d9 —$ :=7r /2-e 



tt/4 

,/ 

sin" $ 

l r /4 r&r 



+ b 2 d$ < 



2 l/2+0 







+ b 2 d§ < 



$2/3 



fo 2 d$ := j!(/3), 



sin$ 



is an increasing 



where bi := max [0 ^ /4 ] = (^V^) since #( $ ) := ^ 

function. Since (0, 1/2) 3 j3 h+ is an increasing function on the interval (0, 1/2), then 

we conclude that 

v(n/2) < j'i(l/2) w -0.8013 := Wl . 



Therefore w intersects = 7r/2 at a negative value. 
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(lb) Lower bound for v(tt/2) in the case u > and 9 <G (7r/4,7r/2). 

Next we calculate a lower estimate for v(tt/2). For this we first observe that 



)l-a 



1 1 

+ cos" 6> + sin" 9 



v 2 > 



1 



1 



2 l/2+0 V cos a Q 



where 63 := min ee [ w /4_ ff /2] M^) = 2. Now we consider the interval J := [7r/4,37r/8] and by a 
direct calculation, we have: 



dv 1 



dO ~ 2V2+/3 V cos Q 



1 



+ 2-v 2 {ir/A) 



> »(ar/8) + jf ^ + d$ ;= ia( fl. 



Since: 



u(3tt/8) := vp(3n/8) 
and since both the functions 



21-20 



2^2 

l + \/2, 



2\/2 
\/2- 1 



(0, 1/2) 9 /? -j. ^(3tt/8) and (0, 1/2) 3 /3 h+ f '* - 

Jo 1 



1 / 1 



/2+/3 V $20 



+ 1 - 2/3+ 1 



are strictly increasing it immediately follows that also j'2 is an increasing function and there- 
fore we can conclude that 

v{n/2) > ji(l/2) := w 2 w -1.6267. 

(lc) The continuation of the right branch through double collision at 9 = ir/2 for u < 0. We 
affirm that w remains negative in the interval J. 

In fact, for u < 



< -W2 1 -" 

d6~ 2\ 



1 



cos«6> 



+ 65 



for 65 := h(3ir/8) and ft. as above. Let 6 6 : = <7(tt/8). Integrating we get: 



(41) 



1 r 7 1 

u(3tt/8) < w(tt/2) + - ^ 2172 

1 Z" 1 / 8 1 
^ v{7T/2)+ 2l 2W 

i 



/2+0 ' 



1 



cos" 6> 



+ 6 S 



d0 < 



$2,3 



+ b 5 



d$ < 



1 z^ 78 1 



$2/3 



+ b 5 



rf$ 



0.0902. 



Therefore 



u> 2 < v(tt/2) < v{6) < v(3tt/8) < 0, for 9 E [3tt/8, tt/2], and u < 0. 
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(Id) From the above inequality we get the estimates w\ > v(9) 2 , for all 6 G [37r/8, 7r/2] and u < 0. 
Using this estimate we can write: 



dv _ 1 
d9 ~ 2* 



)l-a 



1 



1 



1 



V 2 > 



1 



1 



COS ( 



a 7 



sin ft cos" 

where a 7 := min ee [ 37r / 8i7r / 2 ] h(6) — h(n/2) — 2. By a direct integration we have: 

pf/2 



w(3tt/8) > v(n/2) + 



3tt/8 

| 



1 



cos" ft 



a 7 



>v{,/2) + \[ J 2^ 



+ 2 



- tof c?<I> > 



>w 2 + 



tt/8 



$" 



u;| ci<E> > 



>w 2 



it/8 



d<£ w -1.0276. 



(le) The continuation of the right branch through double collision at ft = n/2 for u < 0. We 
affirm that w remains negative at the point 9 — 7r/4. 

In fact, for u < 

cfe 1 r, T 1 

< -a/2 1 -" 

d9 ~ 2\] [cos" 

for fo 5 := h{-K/4) and /i as above. Let b 6 := g(7r/4). Integrating we get: 

f 7T/2 X ' 



(42) 



u(tt/4) < v{ir/2) + 
< «(tt/2) + 



tt/4 



2 l/2+/3 ' 



1 

COS" 



r/4 



1 



1 



tt/4 



2 7 2V2+/3 
1 



$2/3 



dft < 



21/2+/3 



$2,3 



b 5 



d<S> 



0.2237. 



(If) To finish the estimate for the continuation of the right branch through double collision, we 
now consider 9 G [0, 7r/8] and u < 0. We observe that in 9 = tt/8, U{9) has a relative 
minimum and therefore for 9 G [0, tt/8] we have U(9) > U(w/8). Now let v := y/2U(n/8)\ 
thus 

and by a direct integration, we have 

/v\ «(^/8) 1 Z^/ 8 
— arcsin - > - / <i# = 7T/16. 

W „(o) 2 7 

Therefore 



f(0) > 'Csin ( 7r/16 + arcsin 



v(w/8) 



> v sin ( 7r/16 + arcsin 



-0.3741 \ 



0.3810. 



By symmetry, the same arguments holds also for the left branch of W u (c_ ) . Now the thesis easily 
follows by the fact that the flow is of gradient-type with restpoints c_ and c+. □ 
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A. 2 Proof of Proposition |3.6 



The proof of this result will be divided into several steps. 

1. There exists o.q G (1,2) such that for each a G (1, cto) the right branch of W u (c-) 
intersects 9 = it/2 (9 = —ir/2, resp.) with v < and then intersects the section 
v = with an angle 9 G (0, 7r/4). 

2. There exists a\ G (1,2) such that for each a G (ai,2) the right branch of W u (c_) 
intersects = 7r/2 (0 = —ir/2, resp.) with v < and then intersects the section 
v = with an angle 9 G (tt/4, 7r/2). 

3. There exists a* G (ao,ai) such that the right branch of W u (c _) intersects the 
right branch of VK"(c + ). 

In order to prove item 1 we shall proceed as in the proof of the previous proposition. 

(la) There exists a G (1,2) such that for each a G (l,ao), v(tt/2) < 0. 

• Upper bound for v(ir/2) in the case u > and # G (tt/4, w/ 2). 



dv 
d9 



1 



1 



1 



cos' 



sin 



- w 2 < 



1 



2 l/2+/3 



1 



1 



1 



cos° 



< 



< 



1 



21/2+0 



cos" 



where 



Thus 



max h(ff) = l + 2 p , 

8g[k/4,-k/2] 



h{6) 



sm 



tt/2 



1 



7T/4 

1 

2V2+/3 
1 



21/2+0 
tt/4 



u(tt/2) < u(tt/4) + 
= u(tt/4) + 
< v(*/4) + Wj 

= —r-y/i + 2/3+ 1 + — 

2M/2 V 91 



fo 2 — * : =7r/2-0 



1 

sin Q $ 



b 9 rf$ < 



7r/4 



1 



b 2 d$ = 

?r/4 



21/2+0 







bi_ 

$2,9 



where 6i := maxro l7r /4] 



sin$ 



\/2) since g($) := 



b 2 d^^fciCS), 



sin$ 



is an increasing 



function. We observe that (1/2, 1) 3 /3 h-> fei(/3) is a continuous increasing function on the 
interval (1/2, 1). However it changes sign in the interval A\ := (1.4, 1.6); in particular if ao 
is any number in the open interval (1, 1.46136], then we can conclude that 

v(tt/2) < fci(l/2) w -0.1659 := u^. 
• Lower bound for v(tt/2) in the case u > and 9 G (7r/4,7r/2). For this we first observe that 

1 



i 



i 



i 



cos Q 9 sin Q i 



21/2+0 V cos Q 



1 
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where fo 3 := mine e [ w /4. ff / 2 ] M^) = 2. Now we consider the interval J := [7r/4,37r/8] and by a 
direct calculation, we have: 



dv 1 

> 



d9 ~ 2V2+/J V cos 



1 



— + 2-^/4) 



„(,r/2) > ,(3^/8) + ^y]-L- e+ 2-^) d9 > 



Since: 



«(3tt/8) := u^(37r/8) = - 
and since both the functions 
(0 



2 l-2/3 



2\/2 \ ^ / 2^2 \ ^ 



l + \/2 / \\/2-l 



1/2) 9 /3 ^ ^(3tt/8) and (0, 1/2) 9/3^ -jI^-L + l - 2 0+i 

are strictly increasing it immediately follows that also j 2 is an increasing function and there- 
fore we can conclude that 

v(tt/2) > j'i(l/2) := w 2 w -1.6267. 

We observe that for a > 1.1 we have v(tt/2) > ji(0.55) := w 3 ~ —1.5285 Therefore v 
intersects 9 = it/ 2 at a negative value. 

(lb) There exists aoo € (l,ao) such that the continuation of the right branch through 
double collision at 9 = tt/2 for u < 0, remains negative in the interval [7r/4, n/2]. 



In fact, for u < 



1 

dt9 - 2" 





1 






+ &5 


cos«6> 







for 65 := /i(497r/100) and h as above. Let f>6 : = ff( 7r /100). Integrating we get: 



cos" 6> 



d0 < 



(43) 



< v(w/2) + \ 

l /-w/100 1 

2i 2V^ 



tt/100 



21/2+/3 y [$2/3 



' 2 /.. :> - 



$2/3 



&5 



< -0.1659 + 



tt/100 1 

2 1 / 2 +P ' 



$2/3 



w 0.0057, 

which holds for a G (1, 1.1]. Thus the thesis follows by choosing a o = 1-1. Therefore 

w 2 < v (tt/2) < v(9) < w(49tt/ 100) < 0, for (9 e [49tt/100, tt/2], and «<0. 
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From the above inequality we get the estimates w\ > v(9) 2 , for all 9 £ [497r/100, n/2] and u < 0. 
Using this estimate we can write: 



dv 1 L, 

= -\ /2 1 -° ! 

dS 2 V 



1 



1 



1 



1 



COS 



a 7 



sin a 6 cos" 6< 

where a 7 := ^0.0^^/100,1^/2] h{6) = h(n/2) = 2. By a direct integration we have: 

1 



w(49tt/100) > u(tt/2) + 



> 



«(t/2) + ^ 
2 So 



tt/2 

49tt/100 
tt/100 



1 



COS«6> 
1 



+ Q7 



- w\ d9 > 



$< 



+ 2 



tt/100 



> W 2 + 

> 1.1804. 



)l-a 



$a 



+ 2 



d$ > 



• Now we prove that the continuation of the right branch through double collision is negative 
at the point 9 — 7r/4. 



In fact, for u < 



dv 1 
dfl - 2" 



1 

cos Q 



for fo 5 := h(n/4) and /i as above. Let b e := g(n/4:). Integrating we get: 

i r /2 i 



v(tt/4) < 



v(n/2) + \ [ 



< v(tt/2) + 
1 



/4 

7T/4 



2V2+/3 ' 







(44) 



<Wi + 



2 

7T/ ^ 



2 l/2+/3 ' 



cos a 6> 

$2/3 



2 7 2 1 /2+/3 

7T/4 



$2/3 



+ *> 5 



i r 7 i 

^ fei ^ + 2 y 2W 



$20 



< -0.6857 + 



< -0.0713. 



tt/4 



1 



2V2+/3 



$2/3 



+ b 5 



To finish the estimate for the continuation of the right branch through double collision, we 
now consider 9 £ [0,7r/8] and u < 0. We observe that in 9 = n/8, U{9) has a relative minimum 
and therefore for 9 £ [0,tt/8] we have U{9) > U{ir/8). Now let v := ^2J7(tt/8); thus 



dv ^ 1 2 

^2^ V 



and by a direct integration, we have 

sin - 
V v 



arcsm 



w(ir/8) I f*/ s 

> - d9 = tt/16. 
2 Jo 



„(0) 
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Therefore 

v(0) > 5 sin ( 7r/16 + arcsin 



v{tt/8) 



> v sin ( 7r/16 + arcsin 



-0.3741 \ 

5 / 



0.4948. 



By symmetry, the same arguments holds also for the left branch of W u (c_ ) . Now the thesis easily 
follows by the fact that the flow is of gradient-type with restpoints c and c+. 

(lc) There exists aoi £ (aoo, ao) such that the continuation of the right branch through 
double collision at 9 = tt/2 for u < 0, changes sign in the interval [n/A, ir/2] for any 
a £ [aoi,ao). 



In fact, using this estimate we can write: 



dv 1 L, 

= ~\ /2 1 -° ! 

dd 2 V 



1 1 



- V 2 > 



1 



/ 2 l-a 



1 



cos 



+ 07 



sin" 6> cos a 

where a-j := min ee [ 4g7r / 10 o i7r / 2 ] = h(n/2) — 2. By a direct integration we have: 

1 



Wo 



w(49tt/100) > v(tt/2) + 



tt/2 

49tt/100 
tt/100 



1 



1 /-7T/1UU 

>v(,/2) + -J o ^21- 
2 I 



cos c 



1 



+ Q7 



-w\d9> 



ptt/IOO 



> If 2 + 



+ 2 



- wf <i<I> > 



> 1.1804. 



To finish the estimate for the continuation of the right branch through double collision, we now 
consider 9 £ [37r/8,7r/2] and u < 0. We observe that in 9 = tt/4, U(9) has a relative minimum 
and therefore for 9 £ [tt/4, 49rr/100] we have U{9) > [/(tt/4). Now let v := v / 2U{ir/4); thus 



and by a direct integration, we have 



arcsin 



(1) 



dv ^ 1 F2 2 



t>(497r/100) X /-497T/100 



> 



tt/4 



d0 = 3tt/25. 



Therefore 



u(7r/4) > w sin ^3tt/25 + arcsin ^( 3 V 25 ) ^ ^ > £ sin ^3^/25 + arcsin ^ — ^ ~ 0.0511. 

By symmetry, the same arguments holds also for the left branch of W u (c_). Now the thesis 
easily follows by the fact that the flow is of gradient-type with restpoints c and c + . 

(Id) There exists ao* € (aoo,c*oi) such that for the continuation of the right branch 
through double collision at 9 = tt/2 for u < 0, we have v(ir/4) = 0. 

In order to prove this fact we consider the one parameter family of initial value problems: 



(45) 



dv _ 1 
dO ~ 2' 



2 1-Q 




1 


1 




H 


— u 2 




+ sin" 6> _ 






cos" 





i; 2 :=F tt (6» 



, «(7r/4) = - v /2 1 -«(l + 2^ j + 1 ), 
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which can be written in integral form as follows 



v a (6) =v(ir/4)+ / F a (6,v{6))S, V# e (tt/4, tt/2). 

Jtt/4 

If v is a classical ( < ^' 1 )-solution of the ivp above, then it is a "if "-solution of the corresponding 
Volterra integral equation; moreover the function 

/■tt/2 

7:[l,a ) — >R:a — >■ 7(a) := / F a (0,v(6)) d6, 

is continuous, by taking into account the theorem on integrals depending on parameters. Thus 
the function 

V : (1, a ) ->• R : a h-> V(a) := w q (tt/2) := w(tt/4) + 7(a) 

is continuous. Moreover by the monotonicity of the integral of a nonnegative function with respect 
to the domain of integration and by the fact that the map a i-> F a (6,v) increases on (1,2), it 
follows that the map V also increases. Furthermore V(aoo) < and V(aoi) > 0; thus by taking 
into account the theorem of zeros of a continuous and increasing function, there exists a unique 
a* in between such that V(a *) = 0, which means nothing but that v a {n/A) = 0. 
In order to prove item 2, we proceed as follows. 

(2) There exists a\ G (1,2) such that for each a e («i,2), v(ir/2) > 0. 
• Lower bound for v(ir/2) in the case u > and 9 € (7r/4, 7t/2). 
Next we calculate a lower estimate for v(tt/2). For this we first observe that 



1 + 



1 1 

+ 



cos" 9 sin" 9 



,,2 > 



1 



1 



cos«6> 



— V 2 , 



where b 3 := min ee [ T / 4 i7I y 2 ] h{9) — 2. Now we consider the interval J := [37r/8, n/2] and by a direct 
calculation, we have: 



dv 1 
d0-2" 



>1-Q 



1 

COS" 
tt/2 



W 2 (7T/4) 



V (7r/2) > v(3tt/8) + f -± 

/■tt/8 1 



cos° 



2-2«- 1 w 2 (7r/4) d(9> 



1 

2V2+/3 V $2^3 



+ 1-2/3+ 1 d$:=fc 2 (/3). 



We observe that: 

v(3ir/8) := v p (3n/8) = - 
Moreover, the function below 



\ 



2 l-2/3 



1 + 



2s/2 
1 + V2, 



2V2 
V2- 1 



/■Tr/8 -, 

(1/2, 1)3(3^ v p (3n/8) and (0, 1/2) 9 /3 ^ J ^ 



1 

2+/3 V $2£ 



+ 1 - 2/3+ 1 



as well as fc 2 are strictly increasing functions. However it changes sign and in particular it is 
positive in the interval A 2 := [1.7, 2). Therefore by choosing any number ot\ G A 2 we can conclude 
that 

v(tt/2) > fe 2 (1.7) := w 2 w 0.1055 > 0. 

By symmetry, the same arguments holds also for the left branch of W u (c_). Now the thesis 
easily follows by the fact that the flow is of gradient-type with restpoints c and c + . 
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(3) There exists a unique a* G (ao,Q!i) such that v a (n/2) = 0. 
In order to prove this fact we consider the one parameter family of initial value problems: 

i o 



(46) 



dv 



1 



1 



1 



cos" (9 sin Q 6> 



F a (6,v) 



s «(tt/4) = - v /2i-«(l + 2^+ 1 ) ! 
which can be written in integral form as follows 



v a (0) = «(tt/4) + / F a (0,v{0))M, 

Jtt/A 



V6> G (7r/4,7r/2). 



If u is a classical ("if ^-solution of the ivp above, then it is a ^-solution of the corresponding 
Volterra integral equation; moreover the function 



I: [1,2) 



1(a) 



w/2 



F a (6,v(6))dd, 



tt/4 



is continuous, by taking into account the theorem on integrals depending on parameters. Thus 
the function 

V : (1, 2) -> R : a i-» V(a) := v a (ir/2) := v(tt/4) + 1(a) 

is continuous. Moreover by the monotonicity of the integral of a nonnegative function with respect 
to the domain of integration and by the fact that the map a i— > l^ a (0,u) increases on (1,2), it 
follows that the map V also increases. Furthermore V(o.q) < and V(a.\) > 0; thus by taking 
into account the theorem of zeros of a continuous and increasing function, there exists a unique 
a* in between such that V(a*) — 0, which means nothing but that v a (ir/2) =0. □ 



B Some pointwise estimates for the tetrahedral case 



B.l Proof of Proposition 3.7 



The proof of this result will be divided into two main steps. The first correspond to the Newtonian 
case a = 1 and the second to the case a G (0, 1). 

Moreover we refer to the right (left) unstable branch ofW u (e^ 1 ) as that having u > (u < 0) 
in a small neighborhood of in 2? . 

Using uniqueness of solutions, the symmetries and the gradient-like character of the flow with 
respect to v, it is sufficient to prove the following fact. 

(1) The right branch of W u (e^ 1 ) intersects ip — ir/2 (ip — —ir/2, resp.) with v < and 
then intersects the section v = with an angle ip > (resp. ip < 0). 

Proof in the Newtonian case. 

The proof of this claim will be divided into some parts: 

(la) Upper bound for v(tt/2) in the case u > and <p G (<pi, 7r/2), where <p\ = arctan(l/V2)- 



dv 1/1 2V2 

= — A / 1 — V ^ 

dO 2 y cos ip y/i + S i n 2 ~ 



2 V cos <p 
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where 

di := max m(ip) = m {v) '■ 



¥>e[yi,7r/2] ' + sin 2 ip 

We observe that the maximum is achieved at the point ipi since the function m is strictly 
decreasing on that interval; in fact 

,, . 2v / 2sin(/?cos< / 5 
(1 + sin y))' 3 / 2 

Thus 

1 r 72 f~i 

u(tt/2) < + - / ,/ + ^ dp =* :=7r/2 - ¥ , 

2 y Vl V cos <p 



fn/2-tp! 



-J— + dx d<P < 
sm$ 

/ x 1 r ,2 ~ vi fd 2 7 , i 



77/2— (£>i /~J 

6>/6 + / W-j 2 - +di d$ 



/7T 1 \ $ 

where d 2 := max[ 0i7r / 2 - V i] = — - I — — arctan —= I , since the function n(<j>) := — — — 

2 y z v 2/ sin $ 

is increasing on [0,7r/2 — pi]. Since u(pi) = — and by using the quadrature 

J y/^ + bdx = \Jx(a + bx) + ^= log[2Vby/(x(a + bx) + 2bx + a], 
it follows that 

v(tt/2) < Vl := -0.5727 < 0, 
therefore v intersects 9 — n/2 at a negative value. 

(lb) Lower bound for v(n/2) in the case u > and <p e (<pi,n/2). 

Next we calculate a lower estimate for v(w/2). For this we first observe that 



dv i / i 2\/2 , i r~i ; : 

dp 2 y cos p ^/i + sm 2 ^ 2 V cos ip 

where d 3 := min0 £ [ Vl ff / 2 ] m ( < p) = 2. Now we consider the interval J := [7/167T, 7r/2] C 
[pi, it/2] and by a direct calculation: 



dv i rr " ~ : , i r 72 



, >xa/ + 2-v 2 (pi) =► v(7r/2) > v(7/167r) + - / «/ + 2 - v 2 ^) dp > 

dip 2 V cosp 2 J 7 /i 67r V cosp 

1 W 16 /I 
> -1.4994 := u 2 < 0, 

for $ as above and d^ := 2 — u 2 (pi) = 2 — 3-\/6/2; moreover in the last inequality we used 
the fact that sin$ < $. From (la) &(lb) it follows that v 2 < v(6) < v\, for ip E [p 1 ,ir/2] 
and u > 0. 

(lc) The continuation of the right branch through double collision at 6 = n/2 for u < 0. We 
affirm that v remains negative in the interval J. 
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In fact, for u < 



dv 1 / 1 
dO ~ 2Vcos0 



for a 5 := m(7n/16) w 2.0193. Let <i 6 := n(7r/16) w 1.0064. Integrating we get: 

1 r /2 n 

w(3/8tt) < «(7r/2) + - / 4 / +d 5 d<p< 

2 J7n/16 V cos<^ 

( 47 ) , i r /16 /* T „ i r /16 



< wi +0.4723 < -0.1004 < 

Therefore 

«2 < f(7r/2) < v(6) < v(7tt/16) < 0, for e [7tt/16, tt/2], and u < 0. 
Therefore 

v 2 < v{tt/2) < v(6) < v(7n/16) < 0, for (9 e [7tt/16, tt/2], and u < 0. 

(Id) From the above inequality we get the estimates v\ > v(9) 2 , for all 6 € [77r/16,7r/2] and 
u < 0. Using this estimate we can write: 



= Jv/2t%)-^ > 1 J J- + d 7 v\ 
dip 2 2 V cos ip 

where a*? := min ¥ , e [ 77I ./ 16 7r /2] m {f) = f( n /2) — 2. By a direct integration we have: 
«(7tt/16) > V (tt/2) + \ f ' J-L-+d T -vj d0 > 

2 J77T/16 V COS^ 
t/16 



1 P 710 /l 

>v(n/2) + -j^ ) J-+2-vZd* 

i W 16 Ai 

>v 2 + - \ — + d & d$ > t> 3 := -1.0600, 



2 ./n V $ 



where ds : = 2 — f|. 



(le) To finish the estimate for the continuation of the right branch through double collision, we 
now consider ip e [0, 77r/16] and u < 0. We observe that in <p = ipi, U{9) has a relative 
minimum and therefore for ip e [0, 77r/16] we have U(ip) > U(ipi). Now let v := y/2U{ip\); 
thus 

dv 1 

— > 
dip ~ 2 

and by a direct integration, we have 



— — > -y v — v 



in I ^ I > x / d0 = 7tt/32. 

\vJ „(o) 2J 



— arcsm i 
Therefore 

v(0) > usin(77r/32 + arcsin(u(77r/16)/i))) > vsin (7tt/32 + arcsin(u 3 /w)) > 0.1937. 

By symmetry, the same arguments holds also for the left branch of W u (c_ ) . Now the thesis easily 
follows by the fact that the flow is of gradient-type with restpoints c_ and c+. 
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(2) The left branch of ^(ejj intersects cp = —tt/2 with at a negative value of v. 

In order to prove this claim we follow the idea used by |DV99j . 

Now we consider the left branch of W u {e^ 1 ). It is clear that for the same value of tp, v is greater 
along the left branch of W u {e^[ l ) that the corresponding value for the left branch of W u (e^ 2 ), then 
due to the symmetries of the problem it follows that the left branch of W u {e^ l ) meets the section 
v = for some value of p> which is negative. It remains to prove that it actually meets ip = —tt/2 
at a negative value of v. For u < 0, we have: 



dv / . . v 2 

where a(ip) := UTct ^ = — 1 1 ^ If J : = [<^ n _ 1 ,^ n ] and a n := max ve/ii a(ip), 

2 4cos^ 2^/1 + sin 2 

then for tp £ I n , we have: 



By a direct integration, we can conclude that 

(48) v(tp n ) < 2^a^sin [ \\<p n - <f n -i\ + arcsin ' (v " 



rtirn r /t— -l i ■ i /— — 

It is readily verified that a n = for any interval /„ C [— 7r/2, — tp] U [0, (pi\ := A, since ip M> 

a(y>) is positive and non-increasing on that subset; moreover a n = a(<p n ^i) in [7r/2, 7r/2] \ A := B 
since on this set this function is positive and non-decreasing. We define 

(12-nVi 
tp n — — , for each n g 36. 

Solving recursively the inequality, we obtain 

= «(-5v>i/3) < -1.1452. 
On the other side, for p € [—tt/2, 5^i/3], 

dv i r\ ~ 

Ci(/3 2 V COS 

where dg := m(—5ipi/3). Again integrating and letting dio := u{tt /2 — 5<^i/3), we obtain: 

1 r— 5vi/3 / j 
-v(5ipi/a) + v(-ir/2)<~ J +d 9 dp 

2 J- T /2 V cos ip 

(49) ^2 I VshT¥ + d9dCi) 



2 7 V <T> 

< 0.8803 



Now suppose that v{— tt/2) > 0; then from the inequality (49) it follows: 

u(-5<^i/3) > -0.8803 
which is a contradiction. We conclude that v(—ir/2) < for the left branch of W u (e^ 1 ). 



□ 
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Proof in the case a e (0, 1) 

The proof of this claim will be divided into some parts: 

(la) Upper bound for v(ir/2) in the case u > and <p £ (<pi,n/2), where <p>\ = arctan(l/\/2). 



dv 1 1 2 1 -"/ 2 o 1 / 1 2 1 -"/ 2 

- v 2 < 



2 V (2 cos ip) a (l + sinV) Q/2 ~2l/(2cos(^) Q (1 + sin 2 p) Q / 2 



1 2 1+ P 1 / 1 



2 l/2+/3 y cos ^2,3 + (1 + gin 2 ^ ^ 2 l/2+/3 y cog + e l 

where 

2i+/3 

ei := max m(/3, - 2 1 "^ m(/3, := , 2 

v e[ V i,7r/2] (1 + sm 

In fact for each /3 e (0,1/2), the map ip >-> m(f3,ip) is strictly decreasing on the interval 
[<^i,7r/2]; thus the maximum is achieved at the point ip = <p>\. Moreover 

dm 2 1+/3 r, „ , „ .9m 5m 2 2+ ^ sin cj cos p 

¥ = (IW7 [1 ° g2 " l0g(1 + SI11 ^ ^ = ' (l + sin 2 ^ +1 

Thus 

i r /2 / i 

v(n/2) < vfa) + ^j- 2 y] + ei dp =* :=7r/2 - v 

^ rw/2-tpi I ^ 



f 7r/2-^i 



i W 2 -^ 

where e 2 := max[o !7r /2- Vl ] n((3, $) where n(/3, $) = I — — — J . We observe that 



dn _ ( $ \ 2/3 $ 9n _ 2/3 sin $ 

9/3 ~~ Vsin¥y ° S sin $ ' 9$ ~~ $ 



$ \ 2/3 / 1 $cos$ 



sin <f> 7 V sm $ sin $ 



Now since the function [0,7r/2 — </?ll 9 <f> >->• o($) := increases and its minimum is 

sm<P 

achieved at the point $ — 0, then it follows that on that interval is greater than 1. Therefore 
for each $ £ [0, 7r/2— ip\], f3 i-> n(/3, $) is an increasing function. By elementary consideration 
for each /3 e (0, 1/2) the function $ i-> n(/3, <&) is an increasing function. Therefore 



e 2 := 



V6 /?r 1 

arctan 

2 



-i 2/3 ^ 

VE fir 1 

< — arctan — = 

~ 2 \2 y/2 



Now, the function (0, 1/2) 3/5^ Ji(/3) is a negative and increasing function. Therefore the 
supremum is at the point (3 = 1/2. Thus by elementary calculations, it follows that 

v(ir/2) < V! := -0.5727 < 0. 

Therefore v intersects ip — tt/2 at a negative value. 



36 



(lb) Lower bound for v(n/2) in the case u > and p e (ipi,n/2). 

Next we calculate a lower estimate for v(ir/2). For this we first observe that 



dv _ 1 2 
dp~~2V2?P 



1 2 1 +' 3 

+ 



cos 2 / 3 95 (1 + sin 2 Lp)P 



v 2 > 



1 — 

2 V 2 2 ^ 



1 



cos 2 / 3 9? 



+ e 3 



where e 3 := min^g^ w / 2 ] V 3 ) = 2- Now we consider the interval J := [7/167T, 7r/2] C 
[<£i,7r/2] and by a direct calculation: 



efc 1/2 
d^-2V¥0 



cos 2 ? (tp) 
(tt/2) > «(7/16tt) + ^ 



tt/2 

7/167T 



« 2 (</?i) 

22/3 



cos 2 / 3 9? 



/■tt/16 1 /.tt/16 



1 



1 



2/3+1/2 y $2/3 



+ e 4 d$, 



for e 4 := 2 — u 2 ((^i)2 2/3 /2 = 2 — 3^ +1 2 ^; moreover in the last inequality we used the fact 
that sin$ < $. The function 



,tt/16 x / x 

is positive, increasing on the interval [0, 1/2]; moreover j 2 (0) = 0. Thus 

v{ir/2) > -V6 3/ 3 2-3/9 > -1.4993 := v 2 . 
From (la) &(lb) it follows that v 2 < v((p) < V\, for ip £ [ipi,n/2] and u > 0. 

(lc) The continuation of the right branch through double collision at <p = ir/2 for u < 0. We 
affirm that v remains negative in the interval J. 

In fact, for u < 



1/2 



1 



cos 2 ' 3 95 



e 5 



for e 5 := m(/3, 7n/16). Let d 6 := n(/3, 7r/16). Integrating we get: 
(50) 

rn/2 



w(7/16tt) < u(tt/2) + 

< v(tt/2) + 

< -0.5727- 



1 



2/3+1/2 
1 

2/3+V2 



1 



2^/2 / 



Ttt/16 V COS 2 / 3 ip 
/16 

tt/16 



+ e 5 dip < 



J /-tt/16 



|+e 5 d$ 



e 5 d$ < -0.1004 < 



Therefore 

v 2 < v(ir/2) < v{9) < v(7n/16) < 0, for 9 G [7tt/16, tt/2], and u < 0. 
Therefore 

^2 < u(tt/2) < u(0) < u(7tt/16) < 0, for 9 e [7tt/16, tt/2], and u < 0. 
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(Id) From the above inequality we get the estimates v 2 > v(9) 2 , for all 9 € [77r/16, 7r/2] and 
u < 0. Using this estimate we can write: 



2 l+/3 



cos 2 ^ 3 tp (1 + sin 2 <p)P 



> 



1 /_2_ 

2 V 2^ 



1 



COS 2/3 </9 



+ e 7 



,,2 > 



1 /_2_ 

2 V 2 2 ^ 



COS 2/3 (/? 



-e 7 



where e 7 := min ve [ 77r / 16i7r /2] m(/3, y>) = m(/3,7r/2) = 2. By a direct integration we have: 



w(7tt/16) > v(tt/2) + J / 
2 J 7 



r/2 



7tt/16 



'A 

2 2/3 



1 



COS 2/3 (/? 



+ 2 



t) 2 dy) > 



f.TT/2 

/7tt/16 



/A 

2 2 /3 



COS 2 ^ <p 



> v(n/2) + \ [ 

r' 2 i / i , 



■Do ^</9 > 



/ 7 7r/ 

I-tt/32 



>V 2 + 



7o 2^ 



/ 2 V sin 2 ^ $ 



e 8 d$ 



> ^ 2 + ^ V *V + 68 ^ ^3-1-2078 



where e 8 := 2 — u 2 2 2 ^ 

(le) To finish the estimate for the continuation of the right branch through double collision, we 
now consider tp € [0, 77r/16] and u < 0. We observe that in tp = tpi, U{9) has a relative 
minimum and therefore for tp e [0, 7w/16] we have U(<p) > U(tpi). Now let v :— ^2U(pi); 
thus 

dv 1 

~dp ~ 2 

and by a direct integration, we have 



> - \/ v 2 — v 2 



- arcsin 



,„\ "(157T/32) 1 /-W16 

( - > - / ^ = 77T/32. 

w V (o) 2 y 



Therefore 



u(0) > wsin(77r/32 + arcsin(v(77r/16)/w)) > vsin (7tt/32 + arcsin(v 2 /u)) > 0.4183. 



By symmetry, the same arguments holds also for the left branch of W u (c_ ) . Now the thesis easily 
follows by the fact that the flow is of gradient-type with restpoints c_ and c+. 

(2) The left branch of ^"(e^) intersects ip = — tt/2 with at a negative value of v. 

Now we consider the left branch of W"(e 11 ). It is clear that for the same value of <p, v is greater 
along the left branch of W u (e^ l ) that the corresponding value for the left branch of W u (e^ 2 ), then 
due to the symmetries of the problem it follows that the left branch of W"(e 1 " 1 ) meets the section 
v = for some value of <p which is negative. It remains to prove that it actually meets tp = —tt/2 
at a negative value of v. For u < 0, we have: 

dv r~T . v 2 
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where a(/3, ip) :— ^ Tct (^j — 



21+2/3 



1 2 1 +< 3 



cos 2 / 3 <p (1 + sin 2 



If 7„ := [^n_i,^ n ] and a n (/3) := 



max^g/^ a(/3, then for ip € /„, we have: 

By a direct integration, we can conclude that 
(51) v{(p n ) < 2^a„(/3)sin ( \\(p n - <Pn-i\ + arcsin ' 



It is readily verified that o n (/3) = ^^"^ for any interval /„ C [— 7r/2, — y>] U [0, <^i] := A, since 

<^ i— > a(/3,ip) is positive and non-increasing on that subset; moreover a n (/3) = a(/3,<^ n _i) in 
[7r/2,7r/2] \ A := .B since on this set this function is positive and non-decreasing. We define 

(12-nVi 
tp n = — , for each n e 36. 

Solving recursively the inequality, we obtain 

v(ip 3 2) = v(-5cpt/3) < -1.6699. 
On the other side, for ip e [— tt/2, 5<pi/3>], 



dv 1/1 
~dp ~ 2^+V 2 V cos 2 ^ + 9 

where e£ 9 := m(J3,—5<pi/3). Again integrating (by setting $ := tt/2 + ip) and letting dxo ■- 
n(/3, tt/2 — 5</?i/3), we obtain: 

I r— Sipi/3 / J 

— w(5v>i/3) + w(-tt/2) < -ttttj / W - -gs- - + d g dip 



2^ + V2 7_ 7r/2 Y COS 2 / 3 ip 



^ /.7r/2-5ipi/3 / ^ 
< _„,,,„ / ^ / ™ h dq d$ 



(52) " 2^+V2 7 Q V sin 2 ^ $ 

< 0.8721. 



Now suppose that v(—tt/2) > 0; then from the inequality (52) it follows: 

u(-5<^i/3) > -0.8721 
which is a contradiction. We conclude that v(— tt/2) < for the left branch of W u (e^ 1 ). 



□ 



B.2 Proof of Proposition 3.9 



The proof of this result will be divided into several steps. 

1. There exists ao € (1,2) such that for each a € (l,ao) the right branch of W u (c-) 
intersects 9 = tt/2 (9 = —tt/2, resp.) with v < and then intersects the section 
v = with an angle 9 G (0, tt/A). 
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2. There exists a\ € (1,2) such that for each a € {a\,2) the right branch of W u (c_) 
intersects 9 — n/2 {9 = —it/2, resp.) with v < and then intersects the section 
v = with an angle 9 e (7r/4, 7r/2). 

3. There exists a* e (ao,ai) such that the right branch of W u {cJ) intersects the 
right branch of W u (c + ). 

In order to prove item 1 we shall proceed as in the proof of the previous proposition. 

(la) There exists a € (1,2) such that for each a e (l,a ) ? v(ir/2) < 0. 

• Upper bound for v(tt/2) in the case u > and 9 € (7r/4, 7r/2). 



rfw _ 1 
d9~2" 

< 



COS L 



sm 



7 2 < 



2 l/2+0 



1 



cos° 



sin 



< 



1 



2V2+0 V cos a Q 



+ b 2 



where 



Thus 



6»G[7r/4,7r/2] 



sin 



u(7r/2) < v(tt/4) + 
= u(7r/4) + 



tt/2 

tt/4 
1 



1 



1 



2V2+0 y cos a 5/ 

f 7r/4 



fe 2 d# — $ : =7r/2-0 



1 



2 l/2+/3 ^ V sin « $ 

-/ 4 



6 2 < 



+ b 2 = 



= -2^^1 + 2^ + ^ ^ 4 ^+62 d# := fcl (/3) 



where 61 := max[ w / 4 ] I — 



sin$ 



\/2^ since <?(<&) := ( - — — j is an increasing 



function. We observe that (1/2, 1) 9 /3 i-> fei(/3) is a continuous increasing function on the 
interval (1/2, 1). However it changes sign in the interval A\ := (1.4, 1.6); in particular if ao 
is any number in the open interval (1, 1.46136], then we can conclude that 

v(tt/2) < fci(l/2) w -0.1659 := t»i. 
Lower bound for v(n/2) in the case it > and e (tt/A,w/2). For this we first observe that 

1 



/ 2 l-a 



cos" 



Sill 



-V 2 > 



1 



2i/ 2 +0 V cos Q 9 



+ 6 3 -2«-i w 2, 



where 63 := mine e [ w /4_ ff / 2 ] M^) — 2- Now we consider the interval J := [7r/4,37r/8] and by a 
direct calculation, we have: 



dv 1 / 1 ~ ~~ 

T9^¥J^y^9 +2 - v {n/A) 

a 2 1 / 2 

> v(3tt/8) - 



tt/2 

tt/8 x 



2-w 2 (7t/4) d0> 







21/2+0 V $2,3 



+ 1-20+ 1 d<P := j 2 (/?). 
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Since: 

u(3vt/8) := vp(3ir/8) 
and since both the functions 



2 l-2/3 



l + \/2, 



2\/2 
\/2- 1 



(0, 1/2) 3 (3 ^ ^(3tt/8) and (0, 1/2) 3 p ^ j*'* ^L^^ + 1 - 2^ 

are strictly increasing it immediately follows that also j 2 is an increasing function and there- 
fore we can conclude that 

v(tt/2) > ji(l/2) := w 2 w -1.6267. 

We observe that for a > 1.1 we have v(tt/2) > ji(0.55) := W3 ss —1.5285 Therefore v 
intersects 8 = n/2 at a negative value. 

(lb) There exists aoo £ (l,ao) such that the continuation of the right branch through 
double collision at 9 = tt/2 for u < 0, remains negative in the interval [7r/4, n/2]. 

In fact, for u < 



dv 1 /„, 
d6 ~ 2 \ 



cos" 



+ 65 



for b 5 := /i(49tt/100) and /i as above. Let b e := g(n/100). Integrating we get: 



v(49tt/100) < v(tt/2) 



tt/2 



< v(ir/2) + 



2io 



49tt/100 
tt/100 ^ 



2V2+0 



COS' 



dd < 



(53) 



<w 1 + 2 



2 

tt/100 



2 l/2+0 



$2/3 



li$ < 



1 







2V2+/3 y 

tt/100 



< -0.1659 + 



$2/3 
1 



d$ < 



2V2+/3 



J>6_ 

$2/3 



rf$ 



w 0.0057, 

which holds for a e (1, 1.1]. Thus the thesis follows by choosing a o = 1-1. Therefore 

w 2 < v(tt/2) < v(6) < w(49tt/100) < 0, for 6 £ [49tt/100, tt/2], and u < 0. 

From the above inequality we get the estimates w 2 > v(9) 2 , for all 9 £ [49?t/100, tt/2] and u < 0. 
Using this estimate we can write: 



dv _ 1 
dO ~ 2* 



1 1 

1 + ^77T + 



sin" 9 cos" t 



v 2 > 



■)l-a 



cos 



+ a 7 



Wo 
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where a 7 := min 9e [49^/100,^/2] h(6) = h(n/2) = 2. By a direct integration we have: 

r"/2 / 



u(49tt/100) > v(tt/2) + \ [ 
2 J 4 



> 



> w 2 



49tt/100 

1 r\ 



1 

cos Q 6> 



+ a 7 



w\ de > 



t/100 



2 l-a 



$0 



$< 



+ 2 



- u>2 (2$ > 



d$ > 



> 1.1804. 

• Now we prove that the continuation of the right branch through double collision is negative 
at the point 9 — 7r/4. 



In fact, for u < 



dv 1 

dO ~ 2" 



>1-Q 



1 

COS" 



for b 5 :— h(n/4) and h as above. Let 6 6 := c/(7r/4). Integrating we get: 



v(ir/4) < v(ir/2) + 



it/ 4 
it/ 4 



/ 2 1 
21/2+/3 ' 



COS" 



1 r 7 1 

- w(7r/2)+ 2 y 2^ 



/2+/S ' 



$2/3 



de < 



(54) 



7T/ 



1 



21/2+^ 
r 7r/4 



$2/3 



b 5 



1 r 7 1 

^ k ^ ) + 2j 2W 



$20 



+ b 5 



< -0.6857 + 



< 0.0713. 



/■71 

./o 



2V2+0 



$2/3 



(2$ < 



To finish the estimate for the continuation of the right branch through double collision, we 
now consider 6 e [0,7r/8] and u < 0. We observe that in 6 = ir/8, U(0) has a relative minimum 
and therefore for 6 G [0,tt/8] we have U{6) > U{ir/8). Now let v := y /2U{'K /%); thus 



dv 1 2 



and by a direct integration, we have 



Therefore 

v(0) > vsin ( 7r/16 + arcsin 



• ( v 
— arcsin - 

\v 
v{tt/8) 



v(ir/&) 1 f n / 8 
> 

v(o) 2 



dO = tt/16. 



> v sin ( 7r/16 + arcsin 



-0.3741 \ 
v J 



0.4948. 



By symmetry, the same arguments holds also for the left branch of W u (c_ ) . Now the thesis easily 
follows by the fact that the flow is of gradient-type with restpoints c_ and c+. 
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(lc) There exists aoi € (&oo, a o) such that the continuation of the right branch through 
double collision at 9 = n/2 for u < 0, changes sign in the interval [7r/4, it/2] for any 
a e [aoi,ao). 

In fact, using this estimate we can write: 
1 



dv 
~d9 = 2 ' 



)l-a 



1 



1 



1 



sin" 9 cos" 61 



^2\r 



cos 8 



a 7 



where a-j := mine £ [4g 7r /ioo, 7r /2] M^) = h(n/2) = 2. By a direct integration we have: 



w(49tt/100) > v(n/2) + \ [ 
2 J & 

2 I 



tt/2 

49tt/100 
tt/100 



1 



cos«6> 



+ CI7 



d6» > 



+ 2 



tt/100 



> w 2 



> 1.1804. 



To finish the estimate for the continuation of the right branch through double collision, we now 
consider 9 € [3tt/8,tt/2] and u < 0. We observe that in 9 = n/4, U(9) has a relative minimum 
and therefore for 9 e [tt/4, 49tt/100] we have [/(6 1 ) > C/(tt/4). Now let v := ^2U(n/4); thus 



and by a direct integration, we have 



arcsm 



(1) 



t>(49Tr/100) X /-497T/100 



t>(7r/4) 



> 



tt/4 



d6> = 3tt/25. 



Therefore 



, - ■ /"o • M3tt/25)\\ . / -1.1804 
u(7r/4j > u sin I 37r/25 + arcsm I : II > u sin I 37T/25 + arcsm I - 



0.0511. 



By symmetry, the same arguments holds also for the left branch of W u (c_). Now the thesis 
easily follows by the fact that the flow is of gradient-type with restpoints c and c + . 

(Id) There exists an* € (a 00 ,a i) such that for the continuation of the right branch 
through double collision at 9 = n/2 for u < 0, we have v(ir/4) = 0. 

In order to prove this fact we consider the one parameter family of initial value problems: 
(55) (P a ) 



dv _ 1 
d9 ~ 2' 





1 1 




2 l-a 


+ cos" 9 ' sin" 9_ 


— V 2 



k «(7r/4) = -v/2i-«(l + 2/3+1), 
which can be written in integral form as follows 



w a (0) = «(tt/4) + / F a (9,v(9))d9, V e (tt/4, tt/2). 

Jtt/4 
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If v is a classical ^-solution of the ivp above, then it is a "^-solution of the corresponding 
Volterra integral equation; moreover the function 



I [l,a ) 



: a i 



1(a) :- 



r/2 



tt/4 



F a (e,v(6))d6, 



is continuous, by taking into account the theorem on integrals depending on parameters. Thus 
the function 

V : (1, a ) ->• R : a i-> := w q (tt/2) := v(tt/4) + 1(a) 

is continuous. Moreover by the monotonicity of the integral of a nonnegative function with respect 
to the domain of integration and by the fact that the map a t— > F a (9,v) increases on (1,2), it 
follows that the map V also increases. Furthermore V(a 00 ) < and V(aoi) > 0; thus by taking 
into account the theorem of zeros of a continuous and increasing function, there exists a unique 
a* in between such that V(a Qsf ) = 0, which means nothing but that v a (-K/A) = 0. 

In order to prove item 2, we proceed as follows. 
(2) There exists a\ € (1,2) such that for each a e (ai,2), v(ir/2) > 0. 
• Lower bound for v(n/2) in the case u > and 9 £ (7r/4, 7r/2). 
Next we calculate a lower estimate for v(tt/2). For this we first observe that 

i r, 



/ 2 l-a 



1 



1 1 

cos« (9 + sin a 



v 2 > 



1 

COS Q 



— v 2 , 



where b 3 :— mingg^/4 _ 7T / 2 ] h(6) — 2. Now we consider the interval J := [37r/8, n/2] and by a direct 
calculation, we have: 



dv 1 
d9-2" 



1 



cos° 

■ tt/2 



v(n/2) > v(3n/8) + f -± 
> v(3n/8) + 

We observe that: 



v 2 (tt/A) 
IT 



r/8 
r/8 



2 l/2+/3 V cos a 



+ 2-2 a - 1 v 2 (n/4) d9> 



1 



21/2+/3 V $2/i 



1 



+ 1-2/3+ 1 d$:=fe 2 (/3). 



v(3ir/8) := u^(3tt/8) 
Moreover, the function below 



1 



2 l-2/3 



2^ 
1 + V2 



2V2 
V2- 1 



(1/2, 1)3/3^ vp(3ir/8) and (0, 1/2) 3 fi ^ t 

Jo 



r/8 



1 



1 



2V2+/3 V $ 2 ' 3 



+ 1 - 2/3+ 1 



as well as k 2 are strictly increasing functions. However it changes sign and in particular it is 
positive in the interval A 2 := [1.7, 2). Therefore by choosing any number a x e A 2 we can conclude 
that 

v(tt/2) > k 2 (1.7) := w 2 w 0.1055 > 0. 

By symmetry, the same arguments holds also for the left branch of W u (c_). Now the thesis 
easily follows by the fact that the flow is of gradient-type with restpoints c and c + . 
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(3) There exists a unique a* e (a ,ai) such that v a (ir/2) = 0. 
In order to prove this fact we consider the one parameter family of initial value problems: 

I dv 



(56) 



d6 2\ 



1 + 



+ 



1 



cos" sin" e 



-v* :=F a (0,v) 



k «(tt/4) = - v /2 1 -«(l + 2^+i), 

which can be written in integral form as follows 

i-e 

t/4 



v a (0) =v(tt/4:)+ [ F a {0,v{0))d0, V9e (n/4,n/2). 

Jtt/4 



If v is a classical ("if ^-solution of the ivp above, then it is a "^-solution of the corresponding 
Volterra integral equation; moreover the function 



J: [1,2) 



1(a) :- 



it/2 



tt/4 



F a (0,v(0))d0, 



is continuous, by taking into account the theorem on integrals depending on parameters. Thus 
the function 

V : (1, 2) -> R : a i-> V(a) := v a (ir/2) := w(tt/4) + 1(a) 

is continuous. Moreover by the monotonicity of the integral of a nonnegative function with respect 
to the domain of integration and by the fact that the map a i-> F a (0,v) increases on (1,2), it 
follows that the map V also increases. Furthermore V(ao) < and V(ai) > 0; thus by taking 
into account the theorem of zeros of a continuous and increasing function, there exists a unique 
a* in between such that V(ot*) = 0, which means nothing but that v a (n/2) = 0. 
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